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ABSTRACT 

We present the first results from a new 250, 350, and 500 /im Galactic Plane survey taken 
with the Balloon-borne Large-Aperture Submillimeter Telescope (BLAST) in 2005. This survey's 
primary goal is to identify and characterize high-mass proto-stellar objects (HMPOs). The region 
studied here covers 4 deg 2 near the open cluster NGC 6823 in the constellation Vulpecula {I = 
59°). We find 60 compact sources (< 60" diameter) detected simultaneously in all three bands. 
Their spectral energy distributions (SEDs) are constrained through BLAST, IRAS, Spitzer MIPS, 
and MSX photometry, with inferred dust temperatures spanning ~ 12-40 K assuming a dust 
emissivity index [3 = 1.5. The luminosity-to-mass ratio, a distance-independent quantity, spans 
~ 0. 2-130 L Q Mg . Distances are estimated from coincident 13 CO(l — * 0) velocities combined 
with a variety of other velocity and morphological data in the literature. In total, 49 sources 
are associated with a molecular cloud complex encompassing NGC 6823 (distance ~ 2.3 kpc), 
10 objects with the Perseus Arm (~ 8.5 kpc) and one object is probably in the outer Galaxy 
(~ 14 kpc). Near NGC 6823, the inferred luminosities and masses of BLAST sources span ~ 40- 
10 4 L Q , and - 1 5-700 M Q , respectively. The mass spectrum is compatible with molecular gas 
masses in other high-mass star forming regions. Several luminous sources appear to be Ultra 
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Compact H II regions powered by early B stars. However, many of the objects are cool, massive 
gravitationally-bound clumps with no obvious internal radiation from a protostar, and hence 
excellent HMPO candidates. 

Subject headings: submillimeter — stars: formation — ISM: clouds — balloons 



1. INTRODUCTION 

The Balloon-borne Large Aperture Submillime- 
ter Telescope (BLAST) is a 2-m stratospheric bal- 
loon telescope that observes simultaneously at 
250, 350, and 500 usi ng bolometric imaging ar- 
rays |p2caleieEal]2003) . During the first BLAST 
science flight, a 100 hr Arctic flight from Swe- 
den to Canada in June 2005 (BLAST05), BLAST 
conducted the first sensitive large-scale Galactic 
Plane surveys at these wavelengths. The focus 
of these surveys is the earliest stages of massive 
star formation. As noted in the recent review by 
Zinnecker fc Yorkd (|2007l ) our understanding of 
this evolutionary phase is still limited. Over the 
past 15 years submillimeter observations at longer 
wavelengths (350-1200 /zm), e.g. with SCUBA 
on the 15-m Jame s Clerk Maxwell Telescope 
(|Holland et a ll ll999l) or MAMBO o n the IRAM 



199S[ iMotte et al 



30-m telescope (jKrevsa et al. 
20071 ). have opened up studies of the earliest evo- 



lutionary stages of molecular core collapse and 
proto-stellar formation. Such cores lack an inter- 
nal source of radiation and are very cold (< 25 K), 
so that they emit the bulk of their radiation at 
submillimeter wavelengths. Since this emission 
is optically thin in the submillimeter band, ob- 
served flux densities are proportional to column 
density and mass. There is a variety of termi- 
nology us ed for these elusive early stages. For 
example, Zinnecker fc Yorkel ( 2007 ) distinguish 
cores (siz e 0.1 pc) embedded within clumps (size 



0.5 pc). IMotte et alJ (|2007l ) refer to molecular 
cloud fragments which could be high-luminosity 
infra-red protostars, infra-red quiet protostars, or 
high-mass pre-stellar cores. Along the lines of 
the latter, here we refer to high-mass proto-stellar 
objects (HMPOs) as compact sources residing in 
dense molecular clouds, that have the potential to 
form (one or more) massive OB stars, having lumi- 
nosities in the range ~ 1O 2 -1O 5 L , but without 
associated radio continuum emission. This lat- 
ter qualification distinguishes them from massive 
young stellar objects, a later stage in which a hot 
star has formed (perhaps still accreting), provid- 
ing the ionizing radiation necessary to form a high 



emission measure Ultra Compact (UC) H II re- 
gion (and hence radio emission) . Note that at the 
distances of high-mass star forming regions sur- 
veying instruments detect quite massive clumps, 
possibly harboring groups of stars or their precur- 
sors, and not resolving the substructure associated 
with individual nodes of collapse; however, once 
star formation is underway the most massive ob- 
jects dominate the luminosity and ionization. 

BLAST is presently unique in its ability to de- 
tect and characterize cold dust emission from a 
range of pre- and proto-stellar sources, constrain- 
ing the temperatures of objects with T < 25 K 
{(3 = 1.5) using its three-band photometry near 
the peak of the spectrum. An earlier, less 
sensitive balloon- borne submi l limete r telescope 
called ProNaOS (jPupac et al.l 1200 ll) had simi- 
lar goals. Previous studies which do not sample 
the spectral peak of t h e thermal emission (e.g. , 



Johnstone etaD l200d: IPierce-Price et alJ 12 000: 



Kirk et al.l2005tlReid fc Wi lson 200 5tlThompson et all 
2005 [ lEnoch et al.l2006tlHill et al.l2006tlSchneider et al 



200a iThompson et al.T l2006l: lYoung et all 1200^ 
lEnoch et al.1 120071 iMoore et all 120071 ) have been 
limited by their relative inability to measure the 
temperature, producing large uncertainties in the 
derived luminosities and masses. Recent surveys 
with Spitzer MIPS can constrain the tem pera- 
tures of warmer objects ( Carev et al.1 l2005h . but 
the youngest and coldest objects are potentially 
not detected even in the long-wavelength Spitzer 
bands. 

BLAST is also very efficient at mapping, and 
so conducted a series of surveys spread across 
the entire portion of the Galactic Plane that was 
available during the flight. These maps encom- 
pass (primarily) known high-mass star forming re- 
gions: a lOdeg 2 map of Cygnus X; a 4deg 2 map 
in Vulpecula (described in detail in this paper); a 
3 deg 2 map in Sagitta; and finally a 6 deg 2 map in 
Aquila towards the Galactic Ri ng Survey molec- 
ular cloud GRSMC 45.60+0.30 (|Rathborne et al 
|2004[ ). The locations and sizes of these maps are 
shown in Figure [T] The Vulpecula map was stud- 
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ied first, as it is the deepest of the four regions, and 
consists of 6.5 hours of data. The analysis of the 
other three fields, and those from the 2006 Antarc- 
tic flight, including a diffraction- limited 50deg 2 
map of the Vela molecular ridge, will be presented 
in forthcoming papers. 



To put these surveys in context, iMotte et al 



(|2007l ) state that "large fast-mapping images with 
MAMBO-2 are currently the best tool to study 
density structure of molecular clouds with high 
spatial resolution." In about 33 hours of ob- 
servations (and 30 more with the less efficient 
MAMBO) they mapped 3deg 2 in Cygnus X; not 
the entire region but targeted areas of high ex- 
tinction. They catalogued 129 compact sources 
(down to a scale of 15") to a peak flux level of 
80mJy (5 a) and proposed an additional 40 some- 
what more extended structures (2'). Larger struc- 
tures (10') were filtered out. In the survey re- 
ported here, we mapped a contiguous (no thresh- 
old for extinction) 4deg 2 in 6.5 hours. Although 
the very largest scales (including the DC level) 
are absent from our maps, our images maintain 
power well beyond 10'; the images reveal inter- 
esting large-scale structures crossing the entire re- 
gion. As reported below, we catalog 60 sources 
to a flux level of 10 Jy at 250 fim. Extrapolat- 
ing this detection limit using a typical spectral 
shape (T = 20 K and (3 = 1.5) this corresponds to 
170mJy at 1.2 mm. We also fi nd fewer sour c es per 
square degree compared to the lMotte et al. ( 2007 ) 
survey of Cygnus X because this region is slightly 
more distant, not as active, and our survey did not 
target areas of high extinction. There are exciting 
large sensitive surveys to look forward to at these 
same wav elengths with the SPIRE instrument for 



Herschel (jGriffin et al.ll2003l) 



The region chosen for this study was centered 
near NGC 6823 in the constellation Vulpecula 
(£ — 59°), and covers approximately 4deg 2 . This 
region of the Galactic Plane is prominent in im- 
ages of thermal dust emission (e.g., IRAS in Fig- 
ure [J), as well as in the radio and the optical. In 
the radio, the identifiable extended H II region ex- 
cited by the massive stars is called Sh2-86 (nom- 
inal position indicated in Figure [2]). The stellar 
HR diagram for NG C 6823 has been examined by 
Massev et ail (|l995l ). They find an age of 5-7 Myr 
for the bulk of the stars. There are several evolved 
OB supergiants of mass ~ 25 M©. The most mas- 



sive star, 07 V ((f)), with mass ~ 40 M Q , appears 
younger (2 Myr) than the rest. In the optical this 
massive star is responsible for the illumination of 
many "elephant trunks" or "pillars" in the eastern 
portion of the H II region. 

Very high-mass stars evolve quickly into su- 
pernovae. We note that a sup e rnova remnant 
(SNR) discovered bv lTavlor etafl (| 19921 ) is within 
20' of the 07 star in projection (< 12 pc lat- 
erally, assuming the same distance), and even 
closer to one B supergiant - — nominal positions 
are indicated in Figure O iMassey et al. (1995) 
find several pre-main sequence stars of mass 5- 
7Mq, which pass through this stage in less than 
2 Myr. Earlier stages of evolution are revealed by 
the bright Infra Red Astronomy Satellite (IRAS) 
sources clustered in this region. At least seven 
of these objects have been observed in molecu- 



lar line studies (|Beuther et al.l 2002; Zh ang et al 



20051: iBeltran et al J 12009 ) 



This paper describes a technique developed for 
deconvolving the out-of-focus beams common to 
all data from BLAST05 (see iTruch et all l2007h 
and for detecting compact objects in the BLAST 
maps. We report a robust list of 60 submillimctcr 
sources with sizes < 60" diameter that are de- 
tected simultaneously in the three BLAST bands 
(H]). Infra-red photometry for BLAST sources 
is obtained from comparisons with IRAS and 
Spitzer MIPS maps, and the Midcourse Space 
Experiment (MSX) point source catalog (SJ3J). 
The BLAST and infra-red data are combined 
to constrain dust temperatures and integrated 
far infra-red (FIR) fluxes using isothermal mod- 
ified blackbody SED fits (Q. Our analysis has 
benefited from a 13 CO(l — > 0) data cube ob- 
tained as a part of a Galactic Plane survey at 
FCRAO (|Brunt fc Heverl l2007f ). Noting coinci- 
dences between compact structures in the 13 CO 
emission, and comparing velocities and morpholo- 
gies with the VLA Galactic Plane Survey, Spitzer 



GLIMPSE (jWhitnev fc GLIMPSE Team! 120051 ) 
and Digitized Sky Survey images, distance esti- 
mates are obtained for all of the sources in SJU Of 
60 BLAST sources, 49 objects are likely associated 
with the molecular clouds surrounding NGC 6823 
at ~ 2.3 kpc, 10 appear to be in the more dis- 
tant Perseus arm (~ 8.5 kpc), and a single object 
is thought to lie in the outer Galaxy (~ 14 kpc). 
These distances are used to calculate source lu- 
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Fig. 1. — Locations of fields for the BLAST05 Galactic Plane surveys. Results from the Vul pecula field are 



discus sed in detail in this paper. The background intensity image is 100 fxm emission from iSchlegel et al 
(|l998h . 



minosities and masses in Sj6l In this section we 
also present an estimate of the mass function for 
sources associated with NGC 6823. In this survey 
high-mass clumps are seen in a range of evolution- 
ary stages: those with embedded UC H II regions; 
luminous objects without substantial radio emis- 
sion (HMPOs); and cold gravitationally-bound 
clumps (or cores) with low luminosity-to-mass ra- 
tio and no evidence of star formation (0). 

2. BLAST OBSERVATIONS TOWARDS 
NGC 6823 

2.1. Observing Strategy 

BLAST05 is a 2-m Cassegrain telescope, whose 
under-illuminated primary mirror is designed to 
produce diffraction-limited beams with FWHM 
40", 58", and 75"at 250, 350, and 500 /im respec- 
tively. The camera consists of three silicon-nitride 
"spid er web" 



bolometer arrays ([Turner et al 
1200 ll ) almost identical to those for SPIRE on Her- 
" 20031) . with 149, 



schel (jGriffin et al J 120031) . with 149, 88, and 43 
detectors at 250, 350, and 500 fim, organized in 
a hexagonal close-packed pattern. Radiation is 
coupled to the bolometers using 2/A spaced con- 
ical feed-horns, so that the 14' x 7' field-of-view 
(simultaneously imaged by all three arrays) is in- 
stantaneously under-sampled. The telescope must 
therefore scan in order to produce fully-sampled 
images. Scanning is also required to modulate 
the bolometer signals to remove low-frequency 
detector drift. A full description o f the BLAST 



described in lTruch et all (|20071 ). 

The optimal detector noise was obtained by 
scanning the telescope in azimuth at O.ldegs -1 
while drifting slowly in elevation such that the 
scan lines are spaced 65" apart. In addition, ob- 
servations of Vulpecula were split between times 
when it was rising and setting, resulting in cross- 
linked scans at angles of ~ 45°. This cross- 
linking greatly reduces large-scale 1// noise in the 
map. Detail s of the BLAST sc anning technique 
are given in iPascale et al.l (|2007l) . 



2.2. Data Reduction 

The raw BLAST05 data are reduced using a 



telesco pe and detectors is given in IPascale et al 
(|2007l ). The 2005 Sweden flight performance is 



com mon pipeline detail ed in IPascale et al.l (J2007) 
and iTruch et all (|2007l ). First, the 100 Hz sam- 
pled bolometer and gondola pointing data are de- 
spiked and the digital filter responses are decon- 
volved. Time-varying bolometer responsivities are 
tracked using an internal calibration lamp. The 
absolute gain of the instrument (including antenna 
efficiency), determined from regular observations 
of Arp 220, is measured with absolute uncertain- 
ties of 8%, 10%, and 12% at 500, 350, and 250 ^m 
respectively. The relative pointing of the telescope 
within a single map is determined to < 5" rms by 
integrating rate gyroscopes with star trackers pro- 
viding an absolute reference. Residual pointing 
offsets between different observations of Vulpec- 
ula are removed by direct rebinning of the data 
into maps, and aligning the peaks of the seven 
brightest compact sources from each pass. Finally, 
maps are made using a new algorithm called Signal 
And Noise Estimation Procedure Including Cor- 
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relations (SANEPIC). This technique has evolved 
from strategies used to produce maps of the Cos- 
mic Microwave Background. A model is first de- 
veloped for the raw bolometer time-stream data 
in which the astronomical signal is produced by 
projecting the estimate of the map into the time 
domain using the known positions of each detector 
over time, and assuming that the noise is Gaus- 
sian, stationary, and correlated in time and be- 
tween detectors. The brightness of each map pixel 
is considered a free parameter in this model. A 
maximum likelihood solution for the map is then 
found using a direct inversion technique. This cal- 
culation includes the full noise power spectra of 
the bolometers, and propagates knowledge of cor- 
relations in the bolometer data to pixels in the 
map. Low-frequency noise (predominantly slowly- 
varying sky emission) is naturally removed by this 
process. In particular, to reconstruct large spa- 
tial scales in the map, this method gives more 
weight to data taken during a single visit, during 
which the bolometer zero-point drift is minimal. 
The DC level of the map is unconstrained, and 
set to by applying a weak high-pass filter to the 
bolometer data before map-making. T he details 



of SANEPIC are given in lPatanchon et all (|2007n . 



2.3. Overview of the Maps 

The BLAST05 maps of Vulpecula are shown in 
Figure [21 The similarity of the three maps is strik- 
ing, both in the diffuse emission and the peaks, 
demonstrating that the bulk of the measurements 
fall on the Rayleigh- Jeans tail of the thermal emis- 
sion. In this case, they are a good measure of 
dust column density in the molecular cloud and 
the mass of the compact sources. 

The far- infrared (60 and 100 /im) emission 
in this part of Vulpecula is dominated by sev- 
eral luminous high-mass star forming regions. 
Three of the s e hav e be en studied recentl y by 
Beltran et all (|2006l ) and IZhang et all (|2005l ). as- 
sociated with the following IRAS sources (for later 
cross-reference we append the BLAST name from 
Table El): 19368+2239 (V03); 19374+2352 (V05); 
and 19388+2357 (V08). A furthe r four regions 
have been studied extensively by iBeuther et alj 
(|2002j ), using CS multi-line, multi-isotopolog 
observations and the 1.2 mm dust continuum: 
19403+2258 (V18); 19410+2336 (V30); 19411+2306 
(V32); and 19413+2332 (V40). Note that these 
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Fig. 2. — Raw BLAST05 maps of Vulpecula pro 



duced with SANEPIC (jPatanchon et alJ |2007J). 



The similarity between the bands is due to the 
filters sampling primarily the Rayleigh- Jeans spec- 
trum of thermal dust emission. Previously iden- 
tified IRAS sources associated with high-mass 
star forming regions are indicated on the 350 /im 
map (see £12. 3p . It is assumed that most of the 
submillimeter emission originates in a molecular 
cloud complex at the distance of the open clus- 
ter NGC 6823 (~ 2.3 kpc away) indicated in the 
500 fim map (the two circles give the cluster core 
size, and maximum extents, respectively). The lo- 
cations of several OB supergiants, an 07 V ((f)) 
star, and a supernova remnant, all associated with 
NGC 6823, are also shown. Sh2-86 is an extended 
H II region excited by these massive stars. 
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Table 1 
vulpecula map sensitivities. 



Band 


FWHP er b 




OS 


(fim) 


(") (") 


(mJy) 


(Jy) 




Raw maps 






250 


207 88.2 


18.0 


0.32 


350 


206 87.8 


9.9 


0.17 


500 


212 90.4 


9.9 


0.18 




Deconvolved maps 




250 


40 17.0 


350.0 


1.17 


350 


50 21.3 


99.0 


0.41 


500 


60 25.5 


122.0 


0.61 



Note.— Full- Width Half-Power (FWHP) 
is the diameter in which 50% of the beam 
power is contained (equivalent to Full- 
Width Half-Maximum for Gaussian beams). 
The standard deviation (size) of a Gaussian 
beam with the equivalent FWHP is given 
by <7b. The rms noise in a single 18" map 
pixel is given by cr p . The rms noise in the 
equivalent FWHP Gaussian beam is given 
by a s . Values are presented both for "Raw 
maps" produced with SANEPIC ( %T3]) and 
for "Deconvolved maps" f ^2.4[) . 
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do not exhaust the list of bright sources, either at 
submillimeter or far-infrared wavelengths. 

In the submillimeter, the brightest high-mass 
star forming region near IRAS 19410+2336 is vis- 
ible in the north-east part of the BLAST maps 
(Figure [2]). The other bright regions are located 
on a roughly "C" -shaped arc readily apparent in 
the IRAS maps. As with IRAS, the BLAST emis- 
sion is in most cases centrally peaked, representing 
a single massive clump at the observed resolution 
of BLAST. There is also diffuse emission from ex- 
tended clouds of dust, as well as blending of (clus- 
tered) point sources convolved with the beam. 

In the 2005 flight, the point spread function 
(PSF) was not diffraction-limited (40", 58", and 
75" at 250, 350, and 500 ^m respectively) as orig- 
inally designed. The resulting full-width half- 
powers (FWHP, equivalent to the full-width half- 
maximum for a Gaussian) in all three bands are 
approximately 3'5 and exhibit complex shapes 
that vary as a functio n of wavelength a nd posi- 
tion in the focal plane (jTruch et alJl2007l ). 

The beam shape is clearly recognizable at many 
locations across the entire field, suggesting the 
presence of point-like objects with angular scales 
smaller than the beam. However, the overlapping 
beam patterns also indicate that the maps are 
highly source-confused (Figure [2]). Even though 
only 3-8% of the beam power is found in the cen- 
tral diffraction-limited peak at each wavelength, 
the signal-to-noise (S/N) of the maps are high 
enough (compare sensitivity for "Raw maps" in 
Table Q] with typical source flux densities in later 
tables) that we are able to deconvolve a significant 
fraction of the beam to recover angular scales close 
to the diffraction limit. We can then use these de- 
convolved maps to search for compact sources. 

2.4. Image Deconvolution 

Deconvolution in astronomy has a long history 
(see the review bv lStarck et"aT1l2002f ). The convo- 
lu tion problem can b e stated, using the notation 
of lStarck etaD (|2002h . as 



I(x,y) = (0*P)(x,y)+N(x,y), 



(1) 



where / is the observed map, O is the true im- 
age, P is the instrument's PSF, N is measurement 
noise, and "*" is the convolution operator. Invok- 
ing the convolution theorem, this expression can 



be written in Fourier space as 

i(u,v) = 6(u,v)P(u,v)+N(u,v), (2) 

where X is the Fourier transform of A. A naive 
solution to the convolution problem is found by 
simple division in Fourier space, 



k, x Hu,v) a / v . N(u,v) 
0{u, v) = — = U(u, v) + 



P{u,v) 



P{u,v) 



(3) 



However, in practice, this direct inversion 
method amplifies the noise at high spatial frequen- 
cies. Many methods have been developed to solve 
Equation [T] iteratively, including least squares, 
maximum likelihood and wavelet-based methods. 
These algorithms typically require high-precision 
knowledge of the PSF, which we do not have for 
the BLAST05 maps. The effective PSFs vary sig- 
nificantly across the maps due to asymmetries in 
the beam-pattern (which itself is well- measured) . 
Over time the orientation of this pattern on the 
sky changes, such that the effective shape of a 
point source depends on the amount of time spent 
observing at different parallactic angles, and in- 
strumental noise variations. Iterative methods for 
deconvolution were attempted, but the variable 
PSFs proved to be too problematic. We con- 
cluded that strong artifacts due to convolution 
are inevitable and we decided to use the direct 
inversion method, even though it is probably non- 
optimal. 

Equation [3] can be rewritten in real space as 
O = I * K, where K, the deconvolution kernel, 
is the inverse transform of P _1 . In order to sup- 
press the amplification of noise at high frequencies 
we effectively re-convolve the map by a Gaussian, 
G, with width approximately equal to the designed 
diffraction limit of the telescope. This can be com- 
bined with the deconvolution formulation by writ- 
ing K = G/P. Here K is also smoothed at large 
spatial frequencies to further suppress noise spikes 
as P goes to zero. Additionally, spikes in K due 
to zero-crossings in P are clipped. 

The PSFs are estimated by two methods. In 
the first synthetic PSF is constructed based 

on measurements of point sources throughout the 
flight. The synthetic PSF attempts to account for 
the range of parallactic angles over which the field 
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was observed by averaging together appropriately- 
weighted beam rotations. The second method uses 
a point source directly from the map. An isolated 
source from the north-west corner of the map is 
used (V05 in Figure [3]), this having relatively low 
surrounding diffuse emission. A low-order polyno- 
mial is fit to the background and subtracted from 
the source to remove the diffuse emission. The 
deconvolutions are performed using both types of 
PSF. The synthetic PSFs provide the best results 
at 350 and 500 ^m, while the isolated source is bet- 
ter for the 250 /jm map, based on the amplitude of 
the ripples seen in the resulting deconvolved maps. 

In order to reduce edge effects caused by the 
Fourier transforms, both the maps, /, and the 
PSFs, P, are apodized and zero-padded prior to 
deconvolution. The map is apodized over a scale 
of 5' along a rectangle bordering the map and the 
PSF is apodized over 2' along a circle of radius 4'. 

Finally, the variance map produced by SANEPIC 
is propagated through the deconvolution filter, 
providing a noise map for use in source-finding 
and fitting. The noise at each wavelength both 
before and after deconvolution are given in Ta- 
bic [T] under "Raw maps" and "Deconvolved maps" 
respectively. 
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2.5. Compact Source Identification 

We search for individual compact sources in 
the deconvolved maps using source-finding anal- 
ysis based on the use of a "compensated PSF", 
also called the Mexican Hat Wavelet (MHW) tech- 
nique (see e.g., Barnard et al. 2004), which iden- 
tifies objects in confused images by subtracting 
a local background. We apply the MHW with a 
characteristic width equal to the nominal resolu- 
tion (the width of G in the deconvolution) , to the 
BLAST deconvolved maps. The peaks of all ~ 6- 
g clumps in each waveband are identified as po- 
tential sources. Due to the excess noise produced 
by the deconvolution process, many false detec- 
tions are found near bright sources. No robust 
method to automatically reject these false peaks 
was found, since the ripples are not uniform across 
the image. We visually reject any peak that ap- 
pears to be associated purely with the noise rip- 
ples. In addition, 11 of the recovered sources that 
lie slightly further away from the brightest objects 
remain clearly affected by these residual ripples. 
Since the ripples are different in each band the 



Fig. 3. — The grayscale image is a resolution en- 
hanced BLAST 350 /im map, with the locations of 
compact sources (Table [5]) indicated: 49 objects 
associated with the clouds surrounding NGC 6823 
at ^2.3 kpc (circles); 10 objects in the Perseus 
Arm at ~8.5kpc (triangles); and a single outer- 
galaxy object at ~14kpc (diamond) — see 
for a discussion of distance estimates. The res- 
olution enhancement was accomplished with di- 
rect Fourier deconvolution of the PSF (see M2.4|) . 
Residual ringing is an artifact of uncertainties in 
the beam shape which varies slightly across the 
map. 
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net effect is to contaminate the observed BLAST 
colors. In total, ~ 50% of the rejected sources 
have formal statistical significances below 9-ct, al- 
though there is a tail of rejected sources extending 
to ~ 25-cr that have relatively higher systematic 
errors since they lie on the largest ripples next to 
the brightest submillimeter objects. The positions 
of BLAST sources derived from the deconvolved 
maps are shown in Figure [3l 

The list of robust detections in each waveband 
are combined and a circular region around each 
peak in the three maps is then simultaneously fit 
with a single Gaussian. In addition nearby sources 
are fit using multiple Gaussians. In each region, 
a fourth order polynomial baseline is fit to the 
surrounding regions and subtracted. The ampli- 
tudes in each waveband, position, and width are 
all free parameters in the fit. The fit is calculated 
using a non-linear least-squares minimization rou- 
tine. Constraints on source position and width are 
included to reduce the possibility of divergence in 
the fitting procedure. Leaving the source width 
as a free parameter in the fit biases the resulting 
flux densities high, but is necessary since the beam 
varies significantly across the field. This bias is in- 
vestigated in the next section. Given the S/N of 
these sources, however, even after deconvolution, 
the flux density uncertainty is dominated by the 
calibration error. The catalog of source positions 
and BLAST flux densities are given in Table O 

2.6. Monte Carlo Simulations 

In order to estimate the effectiveness of our 
compact source extraction process, a series of 
Monte Carlo simulations are performed. Flux den- 
sities, biases, completeness estimates, and posi- 
tional error distributions are investigated by in- 
serting a point source, convolved with the esti- 
mate of the beam, into the actual BLAST maps 
produced with SANEPIC. The deconvolution and 
source extraction processes are applied and the ex- 
tracted flux densities compared with the inputs. 
The source extraction routine is modified slightly 
from that described in £!2.5l in order to avoid manu- 
ally rejecting unreliable sources in residual ripples 
around bright sources — a simple automatic rou- 
tine was implemented that is able to reproduce the 
manual rejection procedure for the original source 
list to an accuracy of ~ 10%. 

In order to extract the flux density bias and 



errors, 500 iterations are performed at each of 
several input flux densities spanning the range of 
source brightnesses found in the map. At each it- 
eration, a 250 /im flux density is chosen and the 
350 and 500 //m flux densities are assigned based 
on the median source colors found in the field. A 
simulated source is considered detected if an ad- 
ditional peak is found compared to the number 
of sources originally found in the map. In some 
cases, the simulated source will by chance land 
on top of a real source at the same position. In 
such cases, the detection of the simulated source 
is only counted if the flux density exceeds that 
of the real source that it obscures. The result- 
ing output flux density distributions are roughly 
Gaussian with a positive tail. The Bayesian 68% 
confidence limits are found and the distributions 
are fit by the Gaussian that passes through these 
points. We find that the measured flux density is 
biased high by a factor of 1-1.2, independent of 
input flux density (Figure . Based on these re- 
sults, we can describe the spread in flux density 
with a 2-component model: the first is indepen- 
dent of flux density, reflecting noise in the map, 
and the other is linear in flux density, due to un- 
certainty in the fitting process. Positional errors 
are calculated using the same simulations. These 
errors are described by circular Gaussians with a 
ranging from ~ 7" at the bright end to ~ 70" at 
the faint end. Completeness is estimated using the 
same set of simulations, except that the flux den- 
sity fits are not performed. At each of the input 
flux densities, 1000 simulations are performed and 
the number of times that the input is detected is 
counted. The results are shown in Figure [5] 

In total, 60 compact sources are detected in the 
Vulpecula maps. Flux densities and their uncer- 
tainties are given in Table [2j Note that the tab- 
ulated values include a color correction based on 
SED fits; see §4.31 for a description of this proce- 
dure. Figure [3] indicates the positions of the de- 
tected sources. 

3. INFRA-RED PHOTOMETRY 

The BLAST data alone cannot accurately con- 
strain the SED corresponding to thermal emis- 
sion in many of the sources from this survey. At 
temperatures T > 25 K (and j3 = 1.5) all of the 
BLAST filters sample the Rayleigh- Jeans tail. To 
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Fig. 4. — The flux density bias, determined from 
500 Monte Carlo simulations, is shown for each 
waveband. The bias is modeled as a constant func- 
tion of input source flux density and the errors are 
modeled with two components (described in 




100 

S250,in [Jy] 



1000 



Fig. 5. — The completeness as a function of 
250 /im flux density determined from 1000 Monte 
Carlo simulations (described in Q2.6$ . 



measure the luminosity and dust mass, additional 
photometry at slightly shorter wavelengths (~50- 
100 /im) is required to identify the emission peak. 
FIR measurements and upper limits are estimated 
for all of the sources using IRAS 60 and 100 /mi 
photometry and the new MIP S 70 /im map from 
MIPS GAL (|Carev et al.l l2005h . Although most 



of the luminosity is produced at submillimeter 
and FIR wavelengths, a non-negligible portion is 
also observed in the mid-infrared (MIR) (A ~10- 
50 /im). This radiation is emitted from hotter 
dust that lies close to the heat source. Protostars 
produce significant amounts of radiation in this 
band, since the material in the vicinity of the ob- 
ject is optically thick. This material is irradiated 
with UV light which is then re-emitted at longer 
wavelengths. To probe the MIR SEDs of BLAST 
sources we used IRAS 12 and 25 /im measurements 
and the MSxB 

catalog at 8, 12, 14, and 21 /im. 

3.1. IRAS 12, 25, 60, and 100 /im Flux Den- 
sities 

A significant fraction of the BLAST sources 
have clear counterparts in the IRAS Point Sourc e 
Catalog version 2.0 (PSC. lHelou fc Walked Il988h . 
We add the BLAST positional uncertainties in 
quadrature with the semi-major axes of the PSC 
error ellipse to determine search radii for each 
source. The BLAST uncertainties are taken from 
the simulations described in $2M but we conser- 
vatively set the error at the bright end to 30". 
Identifications from the IRAS PSC were found for 
23 of the 60 BLAST sources, as summarized in 
Table El 

For sources that lack PSC counterparts or mea- 
surements in any of the IRAS bands, we produce 
measurements or upper limits directly from IRAS 
maps. We use the IRAS Galaxy Atlas (IGA, 
Cao et al.l[l997l at 60 and 100 i im) and the Mid- 
Infra red Galaxy Atlas (MIGA, iKerton fc Martini 
20001 at 12 and 25 /im). These maps are produced 
using the resolution- enhancing algorithm HIRES 
(jAumann et al.lll990f ) . The PSFs vary across these 
maps, showing strong elongation along the scan 
direction. However, these HIRES maps resolve 
sources that are also detected with BLAST (see 
Figure [6]) yet appear confused in IRAS maps that 
have not undergone resolution enhancement (e.g., 



1 http: //www. ipac . caltech. edu/ipac/msx/msx .html 
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the I RIS maps of iMiville-Deschenes fc Lagache 
2005h . 

We use aperture photometry to measure flux 
densities in the IGA and MIGA maps. A circu- 
lar aperture radius of 2'A is used at 100 /im. For 
the remaining bands an elliptical aperture aligned 
with the scan direction is used, with semi-axes 
1.'8 x 1.'2. These apertures are centred over the 
BLAST coordinates and the maps are integrated. 
Baseline pixel values are estimated as the median 
in an annulus between the outer-edge of the pho- 
tometry aperture and a second circle (or ellipse) 
that is larger by a factor of 1.3. This size was cho- 
sen through trial-and-error as a compromise be- 
tween smaller sizes which encompass fewer pixels 
(baseline estimates with larger statistical errors), 
and larger sizes which suffer greater contamina- 
tion from large-scale extended structure or adja- 
cent sources (baseline estimates with larger sys- 
tematic errors). Since many of the sources remain 
confused, the flux density and baseline measure- 
ments are both checked visually. In cases where a 
small amount of emission is detected in the base- 
line annulus, the confused fraction is excised in the 
estimates. Furthermore, if additional sources are 
seen near the edge of the measurement aperture, 
or if it appears to lie on a bright gradient of back- 
ground emission, the measurement is flagged as an 
upper- limit. Finally, if no source is visible in the 
aperture, it is flagged as a non-detection. The un- 
certainty in all of the measurements is estimated 
as the rms in the non-detections. The results of 
this photometry procedure are summarized in Ta- 
ble m 

To verify our technique, we compare these aper- 
ture measurements with flux densities from the 
IRAS PSC and find no significant bias, with scat- 
ters of 12% at 60 (im and 25% at 100 /im. 

3.2. MIPS 70 /im Flux Densities 

After examining the IRAS PSC, and measur- 
ing flux densities in the IGA maps, 32 of the 60 
BLAST sources lacked a clear detection at either 
60 or 100 /im. To constrain the FIR SEDs in 
these cases, we use mos aics of the 70 /im MIPS- 
GAL (jCarev et al. 2005 ) images downloaded from 
the Spitzer public web site using the Spitzer Pride 
software!! The 8 /jm IRAC image (Figure [T7]) is 



2 http: //ssc . spitzer . caltech. edu/propkit/ spot/ 



also produced using this software. 

Similar to the flux densities obtained from 
IRAS maps in ^3.1( aperture photometry with 
baseline correction is used for the 70 /im map. 
The radius of the measurement aperture is 1' and 
the baseline aperture is an annulus with inner 
and outer radii of 1' and It was noticed that 
these measurements yield flux densities systemat- 
ically lower than the 60 /im observations, despite 
the fact that the SED for every source should 
be brighter at 70 /urn, since these wavelengths fall 
on the Wien tail of the thermal emission in all 
cases. This problem arises because saturation and 
non-linearities affect MIPS 70 /im sources with 
flux densities > 20 Jy (A. Noriega-Crespo, pri- 
vate communication) . Fortunately only the fainter 
BLAST sources require 70 /mi photometry to con- 
strain their FIR SEDs, and adopting a 20 Jy cut in 
the catalog of 70 /im measurements yields useful 
data for 10 BLAST sources with no IRAS detec- 
tions. The uncertainty in this photometry was 
estimated as 3 Jy from the rms of observations 
that were visually flagged as non-detections. A 
comparison between MIPS, IRAS, and BLAST is 
shown in Figure [BJ and the MIPS photometry is 
summarized in Tabled 

3.3. 500-8 /im Map Comparison 

Figure [5] shows an expanded view of a 0?3 x 0?3 
region of the maps described in the previous sec- 
tions, encompassing the BLAST sources Vll, V12, 
V13, V15, and V16. It is clear that the relative 
brightnesses of sources vary greatly over the range 
500-8 /im, due to their different intrinsic SEDs. 
None of the BLAST sources in this region of the 
map have counterparts in the IRAS PSC. Source 
confusion is also an issue; for example, although 
V12 is not in the IRAS PSC, a counterpart is seen 
in the higher resolution IGA 60 /im and the MIPS 
70 /im images (and MIGA as well). V16 has a 
faint counterpart at 70 /im. In the same images, 
IRAS PSC sources not seen by BLAST appear: 
19395+2313 between Vll and V13, 19397+2309 
to the NE of V12, and 19399+2312 W of V16; 
the source SW of V15 is not in the IRAS PSC. 
Based on the brightness and color of 19397+2309 
at 60 and 100 /im (23 and 80 Jy respectively), we 
would expect to have detected this source with 
BLAST. Comparison of the 250 /im image with 
the 70 /im image (where the source is resolved into 
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Fig. 6. — Comparison between BLAST05 decon- 
volved maps, IRAS maps produced with HIRES, 
and the Spitzer MIPS (70 ^m) and IRAC (8 ^m) 
maps. Each thumbnail has a size of 0? 3 x 0? 3. Cir- 
cles on the BLAST maps and IRAC map indicate 
source positions with arbitrary diameters 2'5 and 
2' respectively. The circles and ellipses indicated 
on the IRAS and MIPS maps are the apertures 
that were used for photometry. The SED of the 
coldest object, Vll, is shown in Figure [5] 



two), suggests a detection by BLAST, but it did 
not meet our compactness criteria to be tabulated 
as a point source. At even shorter wavelengths, 
additional sources and nebulosities appear, for ex- 
ample 19397+2315 E of Vll. 

3.4. MSX 8, 12, 14, and 21 Flux Densi- 
ties 

MSX flux densities for the BLAST sources are 
obtained from a cross-correlation wit h the MSX 
Point Source Catalog version 2.3 (jEgan et al 
2003b . As with the IRAS PSC (gSZEJ), MSX coun- 



terparts are identified within a search radius of 
the BLAST positions that varies as a function of 
the source brightness. The positional uncertain- 
ties in the MSX catalog are negligible by com- 
parison (rms ~ 4-5"). We find potential MSX 
counterparts for 40 objects in the BLAST cata- 
log. Inspection of the MSX maps shows that some 
of these sources are found in regions with diffuse 
emission, and, given the angular resolution of the 
BLAST catalog, we cannot presently determine 
whether all of the proposed MSX associations are 
BLAST objects, or other protostars in the same 
star forming region at an older evolutionary stage. 
In cases where a single MSX counterpart is iden- 
tified we assume that it is in fact the same source. 
In cases with multiple proposed MSX counter- 
parts, however, the sum of the flux densities is 
taken as an upper limit. Presently no attempt has 
been made to use the shape of the MSX SEDs to 
justify their associations with the BLAST objects. 
All candidates are given in Table [5] 

4. SUBMILLIMETER MIR SEDS 

With submillimeter-MIR photometry, we fit 
the temperature and bolometric flux produced by 
the coldest dust, as well as the bolometric flux pro- 
duced in the MIR. We note that all of our sources 
are likely composed of regions at different tem- 
peratures, typically a warmer core embedded in a 
colder and less dense medium. Our goal is to use a 
simple SED model as an interpolation function for 
the sparsely sampled photometry to estimate the 
total luminosity of each source, and to infer an ap- 
proximate temperature for the dominant emission 
from cold dust. 

We assume optically-thin emission from an 
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isothermal modified blackbody, 







S V = A[ — ) B V {T), 



(4) 



where A is the amplitude of the SED, B V (T) is the 
Planck function, /3 is the dust emissivity index, 
and the emissivity factor (brackets) is normalized 
at a fixed frequency vq. Note that the amplitude 
can be expressed in terms of a total clump mass, 
M c , the dust mass absorption coefficient kq (eval- 
uated at i/ ), and the distance to the object, d, 



A = 



M c k 
~Rd 2 ~' 



(5) 



Since kq refers to a dust mass, the gas-to-dust 
mass ratio, R, is required in the denominator to 
infer total masses. We adopt hq = 10 cm 2 g~\ 
evalua ted at vq = c/250/im, from Hildebrandl 
(|l983h . A factor of 100 is assumed for R. We 
note that the combined mass uncertainty due 
to k an d R is at least as large as a factor 
~ 2 (e.g.. lHildebrandlll983l : lOssenkopf fc Hennind 
1994t iKerton et al.l 120011 ) , depending on assump- 



tions about the environment, age, chemical com- 
position, and shape of the dust grains. Note that 
our adopted value, along with [3 = 1. 5, is equiva- 



lent to the 1.2 mm opacity adopted bv lMotte et al 
(|2007h 



Equation|3]is fit to all of the available photome- 
try from 500-60 /jm using \ 2 optimization (except 
for some cases where the FIR peak is very broad, 
in which case it is on ly fit from 500-100 /im). As 
in iTruch et al.l (|2007t) the band-averaged flux den- 



sity of the model SED is calculated with knowledge 
of the BLAST filter passbands before comparing 
them with measurements. We also include the cor- 
related calibration uncertainties with the statisti- 
cal uncertainties estimated in £12.61 The data co- 
variance matrix, C, is constructed by placing the 
estimated variance for each data point along the 
diagonal. For the BLAST data points, the diago- 
nals are calculated as the quadrature sum of the 
statistical uncertainties with the calibrat ion un- 
certainties listed in Table 1 inlTruch et al. 



(120071 ). 

The off-diagonal cross-correlation terms are esti- 
mated directly from the Pearson correlation coef- 
ficients in the same table. (3 is not well constrained 
and is therefore fixed to j3 — 1.5 (consistent with 
Figurc[9]) so that only A and T are allowed to vary. 



For the 11 cases where the BLAST colors are un- 
reliable ( §2.5| marked in Table [2]), a temperature 
of 20 K has been adopted. 

4.1. Including Upper Limits in \ 2 

Many of the fainter sources have only BLAST 
detections. In order to make the best use of FIR 
upper-limits, "survival analysis" is used to include 
the upper- limits in the calculation of % 2 . Given a 
measurement that is deemed an upper limit, with 
flux density, s\, and uncertainty, a\, the likelihood 
of the model flux density, s, is obtained from the 
integral of the appropriate tail of the likelihood of 
a detection, 



p{s\s h a\) 



1 



' z2/2 dz, (6) 



(S— S\)/(?l 



where we have assumed Gaussian noise. Equa- 
tion [6] is the "surv ival function" for censored data 
(jlsobe et al.l ll986). The behavior of this function 
is easily seen with a few examples. If the model 
flux density is much smaller than the measured 
limit, then the likelihood goes to 1. If the model 
is equal to the measured limit, the likelihood is 
0.5. The likelihood of models more than a few a 
above the limit drops to ~ 0. 

Noting that \ 2 is the negative log-likelihood 
function, and assuming the likelihood of the up- 
per limits is independent of the detections, then 

X 2 = (s ~ Sd)C _1 (s - «d) T - lnp(s|si, <j\). (7) 

Here Sd are the measured detections associated 
with the data covariance matrix, C . Equation [7] 
gives the expression that we minimize in order 
to fit the region of the data consistent with the 
isothermal SED model (> 60 fim). 

4.2. SED Fits 

Uncertainties for the model parameters A and 
T, and non-linearly dependent quantities such as 
the FIR integrated flux, M c , and their distance- 
independent ratio Lfir/M c , are obtained from 
Monte Carlo simulations. Mock data sets are gen- 
erated from realizations of Gaussian noise, includ- 
ing both correlated and uncorrelated errors, as de- 
scribed by the covariance matrix. The \ 2 mini- 
mization process is repeated for each data set, and 
the resulting parameters and dependent quantities 
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placed in histograms. Means and 68% Bayesian 
confidence intervals measured from the relevant 
histograms are given in Table [6] Note that er- 
rors in T have opposite effects on Lfir and M c , 
which exaggerates the uncertainty in Lfir/M c . 

Figure [7] shows an example SED for V30, the 
brightest BLAST object in the sample. The solid 
gray lines indicate the 68% confidence envelope 
of modified blackbodies that fit the BLAST and 
IRAS 100 /im data. The black line consists of 
the best-fit modified blackbody at wavelengths 
> 100 fini, and a series of piecewise-continuous 
power-laws through the IRAS and MSX data 
points at wavelengths < 100 /im. The best fit tem- 
perature to the submillimeter data is 26 K, consid- 
erably smaller than 46 K estimate d from MIR and 
nun data away from the peak (jSridharan et al " 



2002). The luminosity fraction given in column 5 
of Table [6] is the ratio (as a percentage) of the FIR 
integrated flux from the best-fit modified black- 
body (solid gray lines) to the total bolometric flux, 
estimated from the integral of the black line across 
the wavelength range 2-5000 /im. Similarly Fig- 
ure [8] shows the SED of VI 1, the coldest object 
in the sample (T = 12.3 K). As seen in this exam- 
ple, observations across the BLAST wavelength 
range clearly reveal a turnover in the FIR SED 
for the coldest objects. The lack of counterparts 
in the shorter-wavelength bands is emphasized in 
Figure [6] 

The range of modified blackbody models from 
the Monte Carlo simulations are generally very 
tight (e.g. for V30 AT = 1.0 K, and for Vll 
AT = 1.6 K). The narrowness of this range is 
largely due to the prior constraint that the dust 
emissivity index is precisely = 1.5. There is a 
strong degeneracy between the values of and T; 
for example, at larger values of the steepness of 
the Rayleigh- Jeans spectrum increases, requiring 
lower values of T to compensate for that steepness, 
by bringing the FIR peak to longer wavelengths 
and hence closer to the shape of the data reported 
in the BLAST bands. For the two examples dis- 
cussed here, we re-fit the A > 100 /xm region of the 
SED, letting A, T, and vary independently in 
100 Monte Carlo simulations. For V30 the result- 
ing range of temperatures and dust emissivities 
areT = 27.6±2.8K, and = 1.4±0.3 (compared 
with T = 25.9±1.0K when = 1.5). For Vll the 
ranges are T = 9.4 ± 3.4 and = 2.9 ± 1.5 (com- 



pared with T = 12.3 ± 1.6 when = 1.5). Par- 
ticular realizations from these Monte Carlo sim- 
ulations have been selected to illustrate the ap- 
proximate l-a range of T as dashed gray lines in 
Figures [7] and [SJ Many of the BLAST sources 
have poorly constrained values of T and as in 
the example of Vll. We have therefore chosen to 
fix = 1.5 for the remainder of the analysis in 
this paper, which reduces inferred errors in T by 
factors of ~ 2-3. It also emphasizes the range of 
SEDs uncovered by BLAST, since under this con- 
straint T uniquely determines the wavelength of 
the FIR peak. 

In addition to the effect of the prior on 0, the 
treatment of correlated uncertainties between the 
BLAST measurements lead to inferred tempera- 
ture ranges that are perhaps smaller than what 
one might expect given the sizes of the error bars 
in Figures [7] and [8j For example, V30 has pho- 
tometric errors that are only ~ 2% in all three 
bands (Table [2]). However, the error bars shown 
in these figures include approximately ~ 10% cal- 
ibration uncertainties that are nearly 100% corre- 
lated across the three bands. This type of uncer- 
tainty only affects the absolute scale of the SED, 
A, but not the shape which is encoded in T and 
(for which only the independent 2% photomet- 
ric errors are relevant). For this reason the tem- 
perature uncertainty in V30 (under the model as- 
sumptions) is only AT = 1.0 K, demonstrating the 
superior ability of BLAST to constrain the shape 
of the FIR SED of cold sources when systematic 
uncertainties are accounted for appropriately. 

4.3. Color-corrected BLAST Flux Densi- 
ties 

BLAST filters have broad spectr al widths that 



are ~ 30% of the central frequency (jPascale et al 



20071 ). Since colors sampled by the filters are a 
strong function of temperature for cooler objects 
(T < 25 K), a color correction is required to en- 
able direct comparison with SEDs. Once the SED 
(Equation|4]) has been fit to data by minimizing \ 2 
(Equation [7]), we have the choice of either calcu- 
lating different effective wavelengths for each mea- 
surement, or correcting each flux density at fixed 
wavelengths. We choose the latter, and correct 
the flux densities to precisely 2 50, 350, and 500 /xm 
(this procedure is described in lTruch et al.l[2007h . 
These corrected flux densities are given in Tabled 
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Fig. 7.— Example SED for V30, the bright- 
est BLAST object (T = 25.9 K). BLAST error 
bars include correlated calibration uncertainties 
(~10%). Gray lines show the 68% confidence en- 
velope of modified blackbody models from Monte 
Carlo simulations in which the dust cmissivity in- 
dex is fixed at f3 = 1.5. Dashed gray lines are 
particular realizations from a second Monte Carlo 
simulation in which (3 is also left as a free param- 
eter, here indicating the l-er envelope of T ( §4.2|l . 
The black line shows the best-fit modified black- 
body connected continuously at 100 /im to a se- 
ries of piecewise-continuous power-laws between 
the shorter-wavelength data. 
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Fig. 8.— Example SED for VI 1, the coldest 
BLAST source (T = 12.3 K). Symbols have the 
same meaning as in Figure [7l The FIR spectrum 
clearly turns over in the BLAST bands, and there 
are no detections at the IRAS, MIPS or MSX 
wavelengths (see Figure [6]) . The bolometric flux 
for this source can only be crudely estimated from 
the integral of the modified blackbody. 



15 



and a color-color plot for all 60 sources is shown 
in Figure [9] This plot illustrates the variation in 
temperatures probed by the BLAST wavelengths, 
and confirms that the choice of p = 1.5 for the 
SED fits is close to the center of the distribution. 
We have also tested the color-correction method 
using p = 1 and p — 2 and find that the effect 
on the measured flux densities is small compared 
with the photometric uncertainties. 

4.4. Distance-independent Clump Proper- 
ties 

The range of temperatures inferred from the 
SED fits to BLAST data with reliable colors is 
shown in Figure 1101 The survey has convincingly 
detected sources with a range of FIR peak wave- 
lengths, since the temperature span, 10 K < T < 
40 K, is much larger than the individual uncertain- 
ties that are typically ~ 1-2 K (Tabled . It should 
be noted that these latter uncertainties charac- 
terize the range of plausible temperatures under 
the assumption of the simple SED model (Equa- 
tion d} . There is also a systematic dependence 
on P; for example, increasing p from 1.5 to 2.0 
shifts most temperatures lower by about 5 K, and 
derived masses higher by a factor of ~2. Further 
details on correlations are given in M4.21 

As described above, the luminosity-to-mass ra- 
tio can also be calculated from the SED fits with- 
out any knowledge of the distance. For the simple 
SED model adopted, L Fm /M c cx T 4+/3 and so we 
find a broad range in Lfir/M c ~ 0.2-130 L Q Mq 1 
for the objects with T in the range 12 - 40 K (Fig- 
ure ni 



5. DISTANCES 

To deduce luminosities and masses one must 
have an estimate of the distance to each source. 
There is, however, little information in the litera- 
ture regarding distances to known objects in this 
field. Cross-correlation with the IRAS PSC iden- 
tifies 23 of the 60 BLAST sources with known ob- 
jects (Table |3J) ■ While some of these source s have 
been studied exte nsively ( Beuther et al. 2002t 



Beltran et alj|2006l ). few have unambiguous dis- 
tance information. Important clues are provided 
by the morphology and velocity of the interstel- 
lar material and their inter-relationship. We have 
carried out a multi-wavelength assessment corn- 




Fig. 9. — Color-color plot for the compact sources 
detected by BLAST05. We have omitted several 
sources with unreliable colors, identified in Ta- 
ble [2 and two sources with highly unconstrained 
colors due to faint flux densities. Modified black- 
body models with constant /3, equal to 0, 1, 2, and 
3, increasing towards the bottom, and tempera- 
tures ranging from 5 to 200 K are overplotted as 
solid lines. The dashed line represents the same 
model with a constant temperature of 20 K and 
with p ranging from —1 to 3. The \-a statistical 
error bars are shaded so that the more significant 
detections are darker. The flux densities plotted 
have been band-corrected as discussed in SfLj 
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ther p osition nor velocity (see below and lLockman 
19891 Figure 2) is it part of the main Sagittarius 
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Fig. 10. — Histogram of inferred dust tempera- 
tures, assuming a dust emissivity index [3 — 1.5, 
and excluding 11 objects for which the BLAST 
colors are considered unreliable. The entire distri- 
bution shifts downward by ~ 5 K if instead j3 = 2.0 
is adopted (uncertainties are discussed in ^4.2|) . 



paring the BLAST images and point sources wit h 
IRAS (|Cao et al.l [l997t iKerton fc Martini l2000h : 
MSX; IRAC (gH]); the STScI Digitized Sky Sur- 
vejH (DSS); 21cm radio continuum ima ging from 
the V LA Galactic Plane Survey (VGPS. fstil et alj 
2006); an d spectral line imagi ng in 13 CO (1— >0) 
(FOR AO: iBrunt fc Hevedliooi and H I (VGPS). 

We argue below that 49 of the 60 BLAST 
sources in this field are associated with the molec- 
ular cloud complex within which the open cluster 
NGC 6823 has already formed (Figure W). We are 
then able to adopt the photometric distance mea- 



sured for the stars. We use 2.3 kpc (M assev et al 



1995J), though we note that distances in the range 



1.5 to 3.2 kpc have been reported (IHovle et al 



2003!; lGuetterlll993 IPena et al.ll2003l ). In a pro 



gram measurin g parallaxes using methanol masers 
(jMentenl 120071 ) a distance of 2.20 ± 0.01 kpc is 
found for the m asers associated with V30 (IRAS 
19410+2336; see lSzvmczak et al.ll2000l) . about 0?5 
from the cluster center. NGC 6823 is placed in 
a Galactic context in Figures [Tl] and [12] In nei- 



3 http: //archive . stsci . edu/dss/ 



arm. In addition ten of the sources appear to be 
in the more distant Perseus arm (~ 8.5 kpc), and 
a single object (V07) is thought to lie beyond that 
in the outer galaxy (~ 14 kpc). The sources that 
are not associated with the NGC 6823 molecular 
complex are indicated thus in Table [2] 

The IRAS images (see Figure [T]) show that the 
bright sources and diffuse emission represent a dis- 
tinctive enhancement of size ~ 2°, which as we dis- 
cuss below is also the scale of the coherent molec- 
ular complex. However, within the complex there 
is sub-structure, both in space and velocity. The 
angular diameter of the N GC 6823 cluster is ~ 0?5 
( Kharchenko et al. 2005f) . and this is a part of the 
overall complex. The diffuse radio and optical 
emission from the H II region Sh2-86 (LBN 135), 
the portion of the complex that is being ionized, is 
also of this size. Within this region there are ion- 
ization fronts, clearly associated with the ionizing 
stars in NGC 6823, sculpting the parent molecu- 
lar material. A particularly striking example that 
can be seen in the red DSS image is a silhouetted 
"pillar" pointing at the stars, with an ionization 
front at its end and V39 immediately behind it. 
V43 and V49 are other similar examples. 

At 2.3 kpc, ~ 2° corresponds to a physical size 
of 80 pc, so that all of the parts of the complex are 
at essentially the same distance to within a few 
percent (there could be a systematic error in the 
average distance, as noted), with variations indi- 
cated by slightly different velocities or a silhouette 
indicating relative position along the line of sight. 

To introduce the velocity scale, the radio re- 
combinat i on lin e velocity of Sh2-86 measured by 
Lockmanl (| 19891) in a 9' beam is 29.4 ± 2.1 kms" 1 
with FWHM 24.2 ± 2.9 km s" 1 . This velocity 
width, greater than a typical thermal width, and 
typical of widths found in that survey, could re- 
flect accelerated gas motions in the ionized ma- 
terial and/or be intrinsic to the original neutral 
components comprising this star forming complex. 
It is certainly unrelated to differential Galactic ro- 
tation. For a distance of 2 kpc, the gradient in the 
rotation curve (Figure [T2"]) is ~ 125pckm _1 s, and 
so the velocity spread would in that interpretation 
imply a cloud of extreme elongation along the line 
of sight. Note that this velocity and distance com- 
bination does not fall on the mean rotation curve. 
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There are substantial systematic motions in this 
direction, including considerable emission of both 
13 CO and H I beyond the nominal tangent-point 
velocity (~ 34kms _1 ) and so kinematic distances 
from the rotation curve in this range are unreli- 
able. 




o 

[kpc 



Fig. 11. — Galactic spiral arm model of 
Taylor fc Cordesl (| 19931 ). The dashed line indi- 
cates the line of sight towards 60° longitude. The 
star indicates the location of the open cluster 
NGC 6823 at the adopted distance of 2.3 kpc 
( Massev et al.l Il9951 ). We associate 49 of the 



60 BLAST sources with the molecular complex 
within which NGC 6823 has already formed. An- 
other 10 "Perseus arm sources" are approximately 
8.5 kpc distant (triangle) and one object (V07) ap- 
pears to be in the outer Galaxy (plus). The square 
indicates the tangent point. 



5.1. Velocities 

Using the accurate BLAST source positions, we 
have examined 13 CO (1— >0) spectra of Vulpecula, 
having 46" spatial and 1 kms -1 velocity resolution 
(binned from a native resolution of 0.13kms _1 ), 
respectively. Generally a single compact spa- 
tial/spectral coincident feature is identified (see 
the example of source V12 in Figure [13]). We 
therefore fit Gaussian profiles to the velocity com- 
ponent at the position of each source. The lines 
at BLAST positions have a typical FWHM of 
2.5kms- 1 . Ten of the BLAST/IRAS counter- 
parts have previously-measur ed velocities in other 
spectral lines, specifica lly CS ( Beuther et al.|[ 2002: 



Bronfma n et al.l 11996) and NH^ flM olinari et al. 
199fAISridharan et alJ20"ollZinchenko et alJl997t ) 



Our measured velocities for these objects are on 
average within lkms -1 of the published veloci- 
ties, the maximum deviation being 2 kms - . The 
velocities are listed in Table [2] and are presented in 
Figure [14] The broad peak from 21 to 36kms~ 1 
is similar to that seen in the average spectrum of 
this region, which also includes more diffuse gas. 
The 13 CO velocity of V06, at 43 kms" 1 , still falls 
well within this range. 

Several BLAST sources are in molecular clumps 
that clearly affect the radio emission from the ex- 
tended H II region Sh2-86 seen in Figure [TBI For 
example, V32 is a clear local minimum surrounded 
by an ionization front. The pair V22 and V26 
are also at a minimum. Other objects apparently 
shaping or influencing the steepness of the radio 
contours are V47, V49, and V52. All of these 
sources have CO clump velocities in the range 26- 
31kms _1 , like the radio recombination line veloc- 
ity. 

V23 (with V19 close by) is coincident with a 
bright compact H II region and has the weakest 
CO signature of all the BLAST sources. Addi- 
tional velocity information for V23 was obtained 
using H I absorption against the radio continuum 
source using VGPS data. Velocity components are 
identified at —5, 7.5, 21 (very strong), 27, and 
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Fig. 12. — The solid curve shows the expected 
kinematic distance-velocity relati onship derived 
from t he Galactic rotation curve of lBrand fc Blitz 
(|1993l ) projected along i = 60°. At the adopted 
distance of 2.3 kpc for NGC 6823 we show the 
average radio recombination line velocity of Sh2- 
86 as a star. The dark shaded region indicates 
the spread of 13 CO velocities measured for the 
BLAST sources associated with the main cloud 
complex. The uppermost dotted line indicates the 
velocity of V06, which is higher than the main 
spread but still consistent with systematic motions 
around the mean. The light shaded region indi- 
cates the "Perseus arm sources" and the lowest 
line the source V07 in the outer Galaxy. 
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Channel image of 13 CO at a veloc- 
1 for the region shown in Fig. 6. (6) 



Fig. 13.- (a) 
ity of 29kms~ 
Plot of 13 CO spectrum averaged over 9 pixels cen- 
tered on V12 (solid line) and average over whole 
region in panel a (dashed). The bright clump is 
compact (has a high contrast) both spatially and 
in velocity, making it reasonable to identify it with 
the BLAST source (Table 2). 
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Fig. 14. — Histogram of the 13 CO radial velocities 
associated with BLAST sources. The majority of 
these velocities are very similar to the radio re- 
combination hnev^kocityofthe H II region Sh2-86 
measured bv lLockmanl (|1989l ). shown as the Gaus- 
sian (with arbitrary normalization). The dashed 
line indicates the average spectrum of the entire 
13 CO map, including more diffuse gas. The peak 
at lOkrns -1 is the fairly structureless local gas, 
and that at — 5kms~ 1 is the Perseus arm in which 
we identify 10 sources. 
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Fig. 15. — The grayscale shows the VGPS 21-cm 
radio continuum emission, which dominated by a 
structured H II region Sh2-86 and a shell super - 
nova remnant (diameter ~ 15', Taylor et alj l992). 
Overplotted are the BLAST 350 fim contours, 
which reveal the interfaces between thermal dust 
emission (neutral material) and the ionized gas. 
Circles (sources associated with NGC 6823) and 
triangles (Perseus arm sources) identify the com- 
pact BLAST sources. The squares indicate the 
most massive stars in the open cluster NGC 6823. 
Labeled BLAST sources (V22, V26, V32, V47, 
V49 and V52) lie in molecular clumps that clearly 
affect the radio emission from Sh2-86. 



20 



31kms 1 . T he radio recom bination line velocity 
in a 3' beam (|Lockmanlll989h is -2.8 ± l.Skms" 1 
with FWHM 15.2 ± 2.6 kms _1 (the unusually low 
line width suggests an electron temperature less 
than 5000 K). The spatial position of V23 with re- 
spect to the main CO complex and other IRAS 
and BLAST sources suggest it is at the same dis- 
tance, but the velocities at —5 and — 2.8 km s -1 
indicate a much larger distance ~ 9kpc just be- 
yond the solar circle in the Perseus arm. Nearby, 
V20 and V24 have similarly low velocities, appear- 
ing along a clumpy arc of 13 CO emission at —3 
to +2kms" 1 . V17 and V21, which have 13 CO 
components near ~ Okms -1 that are compara- 
ble in strength to weak components in the nomi- 
nal NGC 6823 molecular cloud velocity range, are 
probably in this distant cloud too. We shall call 
these the "Perseus arm sources." 

The only other BLAST source for which we can 
measure the H I absorption directly is V02. It has 
strong sharp absorption at 31 and 35kms _1 and 
weaker components at lOkms -1 , and possibly —6 
and +43kms _1 as well. The 13 CO velocity of this 
clump is 32kms~ 1 , but there is considerable gas 
near ~40kms _1 in the surrounding cloud on this 
side of the map. 

There is a weak radio source that is coincident 
with V07 and a molecular clump at — 54kms~ 1 . 
This source is too faint in the radio for any pos- 
sibility of measuring H I absorption. We use the 
rotation curve (Figure [T^]) to find a kinematic dis- 
tance of 14kpc, making this an intrinsically lumi- 
nous star forming region in the outer Galaxy. 

On the low side of the velocity distribution, 
there are a few other sources with 13 CO com- 
ponents near ~ Okms -1 that are comparable 
in strength to weak components in the nominal 
NGC 6823 molecular cloud velocity range, namely 
V51, V54, V58 and V59. Using C 18 spectra 



(jBrunt fc Heverl 120071 ) . which probe even denser 
gas, the low velocity component is chosen for V51 
and V58 (these join the Perseus arm sources). V54 
and V59 are seen in projection near the lower out- 
skirts of the NGC 6823 cloud, as depicted in the 
zeroth-moment image (Figure IT6"|) . Since there 
is no compelling morphological evidence to relate 
them to the higher-velocity gas, we also consider 
them to be Perseus arm sources. 

Finally, we note a puzzling result. Using the 
H I absorption technique to break the distance 



ambiguity, and Arecibo observ ations (with a 4 ' 
beam) at the Sh2-86 p osition ( Lockmanl 1989h . 
Kuchar fc Banial ( 1994h find absorption in the 
spectrum out to 50kms _1 and so adopt the far 
distance solution 6.3 kpc. Certainly, there is no 
question about the velocity of the H II region and 
hence the associated molecular cloud and BLAST 
sources. But as we have mentioned above, the ion- 
ization fronts seen on the DSS protruding into the 
molecular cloud are clearly caused by the bright 



OB st ars in NCG 6823 identified bv lMassev et al 
(|l995l ). and these have a photometric distance 
much closer to the near solution. We have looked 
at the VGPS H I data (also with continuum, and 
V spatial resolution) across the face of the radio 
nebula. We find that there are large changes in 
the emission and absorption spectrum, with some 
absorption to 50 kms _1 (in line wings) even at the 
fainter western edge, where the continuum emis- 
sion blends into the (non-zero) background. This 
suggests that there is variable distant background 
continuum emission that is making it difficult to 
apply the technique on this extended source. 

5.2. Morphology 

A map of the 13 CO emission integrated over 
the velocity interval 21-36 km s _1 is shown in Fig- 
ure HU Figure [T7] shows the IRAC 8 /im image 
for comparison. These both trace neutral mate- 
rial, but the latter, heavily influenced by PAH 
emission, requires illumination by sub-ionizing ul- 
traviolet radiation (rather than simply tracing 
column density). Also, 8/xm is a short enough 
wavelength for there to be extinction by large 
column densities, producing t he infra-red dark 
cloud (IRDC) ph enomenon ( Egan et al. 1998; 



Simon et al. 2006). A portion of the molecular 



cloud associated with V03 and V04 produces a 
local minimum in the diffuse 8 /im emission, sug- 
gesting it is in the foreground. On a smaller scale, 
V13 and V56 appear to exhibit similar behavior. 

Overlaid on the figures are the BLAST 350 /im 
contours, which reveal excellent correspondence 
between this measure of the neutral column den- 
sity and the other tracers. Also indicated are the 
BLAST poin t sources and the ionizing stars in 



NGC 6823 (Massey et al. 1995) 



A particularly striking filament, which is a 
prominent dust lane in the red DSS image, is seen 
as a ridge in the BLAST, 13 CO, and IRAC maps 
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Fig. 16. — The grayscale shows the 13 CO emission 
integrated over 21 to 36kms~ 1 . Overplotted are 
the BLAST 350 /im contours which demonstrate a 
close correspondence between thermal dust emis- 
sion and molecular gas in this velocity range. Cir- 
cles (sources associated with NGC 6823) and tri- 
angles (Perseus arm sources) identify the compact 
BLAST sources. The squares indicate the most 
massive stars in the open cluster NGC 6823. Also 
indicated is the position of the supernova remnant, 
and the extended H II region Sh2-86. 
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Fig. 17. — The grayscale shows IRAC 8 emis- 
sion, with the same annotations as Figure [TBI We 
note that the labeled BLAST sources (V03, V04, 
V13, and V56) occur at minima in the diffuse 
8 /im emission, possibly indicating that the sub- 
millimctcr objects are (associated with) infra-red 
dark clouds somewhat in the foreground. A promi- 
nent filament which appears as a dust lane in the 
red DSS image is also indicated (probably related 
to the BLAST sources V13, V15, V16, V38, and 
V41). Optical absorption is also associated with 
V47. These features are discussed in §5.21 
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(indicated in Figure [T7|) . running from V30 to the 
south west (parallel to the Galactic Plane) toward 
V03. BLAST sources V13, V15, V38, V41, and 
possibly V16 appear to be related to this optically- 
dark filament. Other optical absorption is asso- 
ciated with BLAST source regions, such as V47 
projected on the H II region. 

An interesting phenomenon seen in the H I data 
is self-absorption (HISA), caused by colder fore- 
ground material; it is most convincingly found us- 
ing the criteri a of narrow lines an d fine-scale spa- 
tial structure ([Gibson et al.ll2000l ) but can also be 



more widespread. HISA appears at 24kms _1 in 
this filamentary structure. Near 33kms _1 there 
is striking filamentary HISA which extends into 
the V22-V26-V28 portion of the molecular com- 
plex, where the lack of H I emission no doubt also 
reflects a real deficit. At 36kms~ 1 there is more 
HISA which traces the diffuse submillimeter emis- 
sion along the ridge ending at V06. This ridge is 
evident in 13 CO at these velocities, but is not so 
clearly seen with IRAC. 

Examination of the sequence of channel maps 
in the 13 CO data cube, or the first moment (aver- 
age velocity) map reveals a velocity gradient from 
the north-east to the south-west (higher to lower 
Galactic longitudes). Along the long filament the 
velocity changes quite systematically, indicating a 
large scale connectivity between the different por- 
tions of the molecular complex. 

Also revealed in these data are large areas of 
low emission. To the west of the stellar cluster 
is a distinctive void, cleared of molecular gas and 
dust. As expected, because the gas has become 
molecular, there is a general deficit of H I across 
this entire region. This can be seen extending to 
adjacent longitudes in th e longitude-velocity dia- 
gram of lStil et all (120061 ). The deficit is exagger- 
ated adjacent to the H II region (near 33kms _1 
where the gas is molecular) . This can appear more 
localized too, e.g., near V03 and in the V08-V09 
cloud at the same velocity as the 13 CO. 

Another voi d coincides w i th a supernova rem- 
nant found by I Taylor et al. (1992) and is seen in 
finer detail in the VGPS radio continuum image 
(Figure IT5]) . The void is seen in H I and in 13 CO 
as well near 22kms~ 1 . The remnant is not notice- 
ably interacting with the surrounding dense ma- 
terial (e.g., no steepening of the 13 CO contours or 
brightening or deformation of the radio shell) , and 



so perhaps the supernova occured in a pre-existing 
cavity, rather than creating it. The diameter of the 
remnant is only 0?2 (8 pc) and so it could be at the 
average distance of this cloud complex (as the ve- 
locity suggests) but not yet have encountered one 
of its clumpy components. 

6. DUST MASSES AND LUMINOSI- 
TIES 

In $5] we determined distances for all of the 
objects: 49 are associated with the open cluster 
NGC 6823 at ~ 2.3 kpc, 10 with the Perseus arm 
at ~ 8.5 kpc, and 1 object is in the outer galaxy 
at ~ 14 kpc (summarized in Table [2]). Here we use 
these distance estimates to convert the integrated 
fluxes from £j4]into bolometric luminosities, L, and 
to derive the mass, M c , for each source. 

The clump mass M c given in column 3 of Ta- 
ble [6] has been obtained from Equation Note 
that although a mass absorption coefficient at 
250 /im appears in this expression, A is derived 
from the fit to the entire submillimeter and FIR 
data, rather than just the single color-corrected 
BLAST 250 /im data point alone. 

The masses of the clumps discovered range over 
- 15-700 M© in NGC 6823, extending to some- 
what higher masses in the Perseus arm. These 
sources are by definition compact in the submil- 
limeter (and they are compact in CO as well), 
with characteristic sizes D < 0.44 pc at 2.3 kpc 
(taking the de-convolved 250 fim map beam size 
of 40" from Table [I] as a limit) . If observed at 
higher resolution many of these "clumps" would 
pr obably qualify as " cores" in the terminology 
of Zinnecker fc YorkU ( 2007 ). The column den- 
sity for a typical 100 Mq clump, assuming a 
uniform sphere limiting diameter D, is of order 
10 23 (0.44pc/£>) 2 cm~ 2 (visual extinction to the 
cloud center Ay ~ 22(0.44pc/D) 2 and density of 
order 10 5 (0.44 pc/D) 3 cm~ 3 ). The Bonner-Ebert 
critical mass for a gravitationally-bound clump, 
above which collapse of an unstable configuration 
occurs, i s given by M rrit /M ^ « (D/pc)(T gas /K) 
(see e.g., Kerton et al.ll2001 . Equation 4). In these 
dense clumps, T gas is probably close to the dust 
temperature and so for the derived T and M c 
most of the clumps that we have identified would 
be unstable and are likely supported by additional 
pressure from turbulence or magnetic fields. On 
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release of this pressure, they have the potential to 
form a massive star or even a small cluster. Inter- 
estingly, several of the lowest mass clumps would 
be just stable if they were as extended as 1 pc. 

The range of luminosities does indeed extend 
to 10 4 Lq in NGC 6823, which corresponds to the 
luminosity of a BO. 5 zero age main sequence star 
(|Panagialll973j ). But as we shall discuss further, 
at the low end, ~ 40 L , there is by contrast no 
indication of significant star formation. 

Luminosity as a function of mass is shown in 
Figures HU and [HI This is essentially a reparame- 
terization of the BLAST flux density vs. tempera- 
ture plane, in terms of intrinsic source parameters. 
Because of the different distances for the two main 
groups, the Perseus arm sources and those associ- 
ated with NGC 6823, the scaling of the axes is 
different, requiring similar but separate diagrams. 
Recall that the derived quantities are based on a 
SED fit assuming a modified blackbody SED with 
[3 = 1.5. Using this same model, the dotted lines 
are loci (with varying T) of constant 250 fim flux 
densities ranging from 10 to 1000 Jy. Orthogo- 
nal to these, the dashed lines are loci (with vary- 
ing observed flux density) at constant T = 10, 
20, 30, and 40 K. With our adopted model, these 
correspond to constant L-pui/M c , values of 0.07, 
3.2, 30, and 140 L Q Mq 1 , respectively. The 11 ob- 
jects with unreliable BLAST colors ( §2.51 marked 
in Table \2§ have been assigned a temperature of 
20 K and so are seen distinctly along that locus. If 
these sources have different temperatures (proba- 
bly lower given the usual lack of IRAS PSC coun- 
terparts) , they would move along the constant flux 
density locus (down and right if cooler). These 
Figures can be used to investigate which early 
stages of star formation might be present. To this 
end, sources have been marked which lack IRAS 
PSC or MSX counterparts, or are UC H II regions 

6.1. Mass Spectrum 

The clumps that we have detected are much 
more massive than stars and therefore represent 
an early stage in star formation that is not easily, 
even empirically, related to the stellar initial mass 
function. Fragmentation and efficiency affect the 
transfer function in unknown ways. Nevertheless, 
it is interesting to investigate the mass spectrum, 
which we have done for the 49 BLAST sources as- 
sociated with the NGC 6823 cloud complex. The 
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Fig. 18. — Luminosity vs. mass for the 49 objects 
believed to be associated with NGC 6823. Un- 
certainties in mass and luminosity are estimated 
from the range of temperatures consistent with 
the BLAST and FIR data under the assumption 
of the simple SED model of Equation |4j Not in- 
cluded are additional sources of error involving 
the distances to the objects and adopted values 
of (3, k , and gas-to-dust mass ratio. The dashed 
lines are loci at fixed temperatures (and L-pm./M c ) 
of 10-40 K, assuming a modified blackbody SED 
with (3 — 1.5. Orthogonal to these are loci (dot- 
ted lines) of constant 250 /im flux density, rang- 
ing from 10 to 1000 Jy, using the same model. 
Sources lacking IRAS PSC or MSX counterparts 
are marked with circles and squares respectively. 
Ultra Compact H II region and Hyper Compact 
H II region candidates are indicated with single 
and double triangles, respectively. 
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Fig. 19. — Same annotations as Figure [T8l but for 
the 10 objects that are associated with the Perseus 
arm. 

masses of individual sources were placed in log- 
arithmically spaced bins (2 per decade); a lower 
limit on the error is then estimated from the Pois- 
son uncertainty for each bin. The resulting mass 
function is shown in Figure I2TJ1 

At masses > 100 M Q , the distribution is consis- 
tent with the mass function s of m olecular clouds. 
For example, Kramer et al. ( 19981 and references 
therein) found a power-law index a ~ 1.7 in 
a mass range which encompasses the mass in- 
terval considered here. For reference, we have 
shown their power-law scaled vertically such that 
it passes through the high-mass bins in Figure l20l 
dN/dM = 150(A//50M Q )- 1 - 7 . Note that adopt- 
ing a different value of dust emissivity affects 
the inferred dust mass, since it is highly corre- 
lated with the fitted temperature ( ^4.2|) . and hence 
mass. For example, choosing (3 = 2.0 decreases 
temperatures by ~ 5 K, and increases cloud masses 
by a factor ~ 2. This effect shifts the mass func- 
tion to the right in this figure, but does not affect 
the power-law index. 

We note that the mass function in our sample 
appears to turn over at masses below ~ 100 M Q . 
However, the completeness of the survey also 
drops at these masses, with uncertainties that 
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Fig. 20. — Mass function of the 49 objects as- 
sociated with NGC 6823. Poisson error bars are 
shown, although there are additional sources of 
uncertainties for the parameters of Equation [5] for 
each object. For reference, the dashed line is a 
power- law, 150(Af/50M o )- 1 - 7 . The two lowest 
bins suggest a turnover, but this shape is prob- 
ably due to incompleteness. To illustrate this 
the 250 jLim flux density completeness (Figure [5|) 
has been converted to a mass completeness us- 
ing Equation [5] assuming various temperatures 
T = 10, 15, and 20 K. These are shown as dotted 
lines, arbitrarily scaled such that 100% complete- 
ness is 40 objects per bin. 
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are difficult to quantify given the range of ob- 
ject temperatures sampled (Figure [TSf . To at- 
tempt to assess the mass completeness, we con- 
verted the x-axis of the flux density completeness 
function (Figure [5]) into mass using Equation [5j 
assuming consistently a distance of 2.3 kpc and 
kq = 10cm 2 g _1 . Some results are shown in Fig- 
ure [20] If we assumed that the sources potentially 
missed had a temperature of 20 K, then complete- 
ness has a negligible effect, but for lower tempera- 
tures, which seem more probable, the effect is sig- 
nificant, because lower temperatures correspond 
to higher masse s for a give n flux densit y. See 



Johnstone etaD (|2000h and lEnoch et all (|2006h 



for further discussions regarding the completeness 
of mass spectrum estimates from submillimeter 
surveys. If, however, the turnover were real at 
~ 100 Mq, our survey may be probing an in- 
teresting effect dis cussed by other authors (see 
Beltran et al. 20061 . and references therein). They 
claim that the point at which the mass distribu- 
tion flattens is proportional to the mass range of 
clumps in star forming clouds. 

7. DISCUSSION 

7.1. Stages of Pre-stellar Evolution 

This particular survey was not unbiased, since 
it targeted a region which has already formed mas- 
sive stars and which contains bright IRAS sources 
which are protostar candidates. But it does fulfill 
one of the main goals of BLAST, to explore the 
sites of formation of high-mass stars and to exam- 
ine the pre-stellar evolutionary sequence. In par- 
ticular, with the BLAST wavebands we are sen- 
sitive to the earliest phases of massive star evo- 
lution, which would be molecular clumps lacking 
an internal energy source and thus cool with low 
Lfir/M c . In Figure fT8l there are indeed many such 
low luminosity-to-mass ratio sources that are too 
faint and cool to have IRAS and MSX counter- 
parts. 

An interesting perspective is that the luminosity- 
to-mass ratio for diffuse cirrus in the local in- 
terstellar radiation field is of order ILqMq . 
The low ratio (and lower temperature) found in 
these molecular clumps, even in the presence of 
a stronger impinging radiation field, is a result of 
the tremendous opacity to this external radiation, 
uncompensated by an internal source. In abso- 



lute terms the luminosity of a cool clump can be 
quite impressive, because of the large mass, and 
yet not be a reliable indicator of there being star 
formation within. 

Once a high-mass pre- or proto-stellar object 
begins to form, there is considerable luminosity 
to be reprocessed, resulting in higher Lfir/M c 
and warmer dust. BLAST can of course see such 
objects as well, even those without MIR associa- 
tions. There is not yet a unique definition or set of 
physical properties associated with HMPOs, but 
there is general agreement that they are heavily 
obscured objects (even making detection difficult 
in the MIR) that have not yet reached the stage at 
which the emergent Lyman continuum is sufficient 
to ionize the surrounding medium to form an H II 
region, and so they are undetected in the radio 
continuum as well. With spectral energy distribu- 
tions peaking at wavelengths > 200 /im they can 
be extremely weak and undetectable by IRAS at 
25 /i m, thus excluding them from previous surveys 



(e.g., Sridharan et al. [20021 ) . Therefore the most 



unique BLAST objects, which are good HMPO 
candidates, are those with no IRAS or MSX coun- 
terparts, and no radio emission. 

In the subsequent phase, when the massive pro- 
tostar has become hotter, near the zero age main 
sequence, the ionizing capability is enhanced and 
a small region in the dense surroundings is ion- 
ized. Historically, such "ultra compact" (UC) H II 
regions have been the primary means to inves- 
tigate the early phases of OB stellar evolution. 
This is a productive approach, since it deals with 
intrinsically luminous objects which are readily 
detectable, and the radio emission provides un- 
equivocal evidence for an ionizing star. The more 
evolved, and therefore warmer, sources are more 
likely to have a MIR counterpart and be detectable 
by MSX and/or IRAC, than the colder and/or 
more deeply embedded sources. We do see such 
objects, as identified in Figures [15] and \T§\ (see 

In the past decade it has become clear that 
"hot-cores" represent a stage in the massive star 
formation process, probably somewhat earlier 
than UC H II regions (see, e.g.. iKurtz et al. 2000, 
and references therein). The precise details of 
the transition are still lacking, and our under- 
standing is complicated by the co-existence of ob- 
jects at different evolutionary phases within the 
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same cluster. Also, the distinction between hot- 
cores and UC H II regions is often dependent on 
the sensitivity of the available observations, as 
several hot-cores are found to emit radio contin- 
uum emission if observed with enough sensitiv- 
ity and angular resolution (opacity can also be 
a factor for dense objects). Likewise, although 
some authors list hot-cores as a subset of HMPOs 
(e.g., ISridharan etaD 12002) . hot-cores are gener- 
ally warmer, with Tk > 100 K, whereas objects 
called HMPOs usuall y have temperatures much 



lower than 100 K (e.g. jMolinari et al.lll996l ). Also 



as noted bv lSridharan et ai* ( 20021 ). among higher 



luminosity objects, a low luminosity-to-mass ratio 
relative to that for U C H II region s migh t indi- 
cate a HMPO, though [Beftran et all (|2006h found 
no evidence to support this difference between 
their "Low" and "High" (less and more evolved) 
sources. 

None of the objects that we detect has an in- 
ferred temperature greater than 40 K (assuming 
(3 — 1.5), even those that are UC H II regions. The 
temperature that we derive from the SED charac- 
terizes the bulk of the dust emission (more than 
50 %) but of course is a simplification. When a 
protostar forms, the dust nearby is hotter than in 
the parent molecular clump, and also appears as 
a more compact source. Thus what is measured 
and interpreted as mass and luminosity depends 
on the radiative transfer and temperature gradi- 
ent as well as the angular resolution of the mea- 
surements. While L reprocessed by dust (and the 
amount of ionization gauged from the radio con- 
tinuum emission) are calorimeters for the embed- 
ded stellar population that is forming, the mass is 
that of the entire parent clump, depleted by ac- 
cretion onto stars and later cloud dispersal. The 
luminosity-to-mass ratio can therefore be differ- 
ent, even for the same stellar population, being 
maximized for optical depth of order unity in a 
compact (and therefore warm) circumstellar re- 
gion. 

A plot like Figure [18] could be used to exam- 
ine the evolutionary sequence, whether at constant 
M c or otherwise, determining the lifetimes for each 
stage from the relative numbers of objects. How- 
ever, this is complicated in practice by effects like 
clustering and detection bias and the fact that star 
formation does not proceed at precisely the same 
rate everywhere in a molecular clump. Thus with- 



out adequate resolution it would be difficult to 
discriminate between material associated with an 
UC H II region and another nearby cold HMPO, 
for example. This problem can be remedied by 
high-resolution and high dynamic range observa- 
tions (as anticipated with Herschel and ALMA); 
indeed in some cases intcrfcrometric observations 
have led to the detection of HMPOs i n the same 
fields as the UC H II regions (see, e.g., Olmi et all 
2003, and references therein). Likewise, objects 



in Figure [18] with candidate IRAS and/or MSX 
counterparts, might in fact harbor a HMPO that is 
associated, but not coincident, with a MIR source 
at a later evolutionary stage in the same cluster. 

7.2. Ultra Compact H II Regions 

We have detected several bona fide UC H II re- 
gions, by which we mean compact objects whose 
bolometric luminosities and radio continuum emis- 
sion are consistent with an embedded massive star 
near the zero age main sequence. Although within 
our beam there might be stars forming over a 
range of masses, the luminosity (and even more 
so the ionizing radiation), is dominated by the 
most massive and hottest stars with the earliest 
spectral types. Therefore, while not precise, it is 
still useful to characterize the region using a sin- 
gle spectral ty pe. We u se luminosities and ioniz- 
ing rates from Panagia ( 19731 ). which to the accu- 
racy required here are close to more modern values 
(jSchaerer k de KoterlH997h . 

These sources can be been identified as UC H II 
candidates using colo r criteria based on IRAS me& - 
Wood k Churchmen] (Il989h . 
Bronfman et al.l (| 19961 ) and 



surements a s in th e 
Kurtz et all (Il994l). 



Watson et al.l ( 20031 ) catalogs (qualitatively, these 
sources are relatively weak at 12 pm). Using the 
BLAST data we are able to find very accurate tem- 
peratures. Radio continuum measurements in the 
survey region at 1.4 GHz are avail able from the 
NRAO VLA Sky Surv ey (NVSffi ICondon et al 
19981) and the VGPS (jStil et allbood) . In ad 
dition, we searched the second and third MIT- 



Gree n Bank 5 GHz "MG" s urveys (|Langston et al 
19901: iGriffith et all Il990h which, though much 



shallower than NVSS, can provide a second flux 
density to estimate the spectral index between 
1.4 and 5 GHz. Some additional targeted obser- 



*http:/ /www. cv.nrao.edu/nvss/ 
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vat ions at 4.9 GHz u sing Arecibo are described in 



Watson et al 



(120031). 

We will discuss these sources briefly in order of 
decreasing distance. V07 is the source in the outer 
Galaxy, which though fairly luminous (3 x 10 4 L© ) , 
is not among the more conspicuous sources in the 
maps. It is consistent with excitation (both lumi- 
nosity and ionization) by a BO star (only 20 M Q 
as compared to 600 M Q for the clump). In the 
VGPS radio image this source is extended and in 
the NVSS it is double. The second component, to 
lower longitude by about 1', coincides with a sec- 
ond IRAS source which appears at 12 and 25 fixsx. 

V23, in the Perseus arm, is the most luminous 
source in this BLAST survey (7 x 10 4 L Q ; Fig- 
ure I19p . and a prominent radio, submillimeter, 
and infra-red source, despite its distance. Its ex- 
citation is consistent with a star as early as 07-8 
(30 M©, as compared to 500 M Q for the clump). 
The emission is extended, overlapping a second 
BLAST source V19, and is no doubt more com- 
plex than this simple description. 

V02 and V05 are among the sources to which 
we have assigned the NGC 6823 distance (Fig- 
ure ITS)) . The luminosity implies a B2 zero age 
main sequence star, but the ionization (50mJy ra- 
dio continuum) requires B0. 5 (for these relatively 
cool stars the ionization is extremely sensitive to 
the temperature). A further note is that the radio 
source, while being aligned with the IRAS source, 
is offset by 30" from V02. 

V30, the brightest BLAST source and the 
most luminous source at this distance, is by con- 
trast not an NVSS/VGPS source at 1.4 GHz. 
However, using the V LA B array at 8.2 GHz, 



Sridharan et al 



(2002) found weak emission at 
the level of 1 mJy. Lowe r resolution Arecibo mea- 

" 20031) 



surements at 4.9 GHz ([Watson et al 
3mJy. Relative to V02 and V05 



give 



V30 is under- 
luminous in the radio continuum by an order of 
magnitude. This could be because the central ob- 
ject, though luminous, has not become hot enough 
on approach to the zero age main sequence (it 
is a HMPO), or because of opacity if the sur- 
roundings of the protostar are very dense. We 
note that about 10% of th e sources observed by 
Wood fc Churchwell ( 19891 ) are optically thick at 
2 cm, and a higher percentage may be applicable 
for the BLAST sources, since they probe the earli- 
est stages of massive star formation. V30 might be 



wha t is called a Hyper Compact (HC) H II region 
(e.g.[Keto 2003), and has been marked accordingly 
in Figure [HI 

The situation for the other comparably-luminous 
sources V08, V18, V32, and V52 is less clear. All 
but V18 meet the color criteria for being an UC 
H II region. However, none are NVSS sources. 
The detection of the V18 and V30 in targeted 
Arecibo observations ( Watson et al.|[2003h appears 
to be contradicted by the lower upper limits of 



1 mJy with the VLA B array (jSridharan et al 



2002). This might suggest contamination of the 



Arecibo beam by diffuse continuum emission; as 
noted in §5.11 V32 is projected against the bright 
patchy Sh2-86 emission. On the other hand, V18 
is coincident with a small (20") cometary opti- 
cal nebula and the Arecibo flux of 8mJy might 
be the more relevant if the VLA resolved the ex- 
tended emission. V18 may be another example 
of a HC H II region, and it is marked as a can- 
didate in Figure [18] For V08 and V52 there are 
no sufficiently sensitive observations to isolate any 
weak emission, or alternatively, to qualify them as 
HMPOs. 

Most of the radio sources in this survey re- 
gion are not in fact associated with BLAST and 
IRAS sources (and vice versa). One must there- 
fore be cautious of chance coincidences (for ref- 
erence the surface density at the ^lmJ y level is 
< 0.007 arcmin" 1 . iFomalont et al.lll99ll ). There 
is a VGPS/NVSS source offset by 75" from V40 
and its I RAS counterpart, al so detected faintly at 



Arecibo ([Watson et al.ll2003[ ). The level of radio 



emission is at least an order of magnitude greater 
than would be e xpected for the relatively low lu- 
minosity B2 star ( Panagiall973 ). so that it is prob- 
ably not as sociated at all. U s ing th e VLA B array 
at 8.2 GHz, Sridharan et alJ ( 20021 ) placed a more 
definitive upper limit of 1 mJy. Note that such 
a HMPO will never become a significant H II re- 
gion unless the stellar mass is increased through 
further accretion. There are NVSS sources near a 
few other even lower luminosity BLAST objects, 
V35, V39, V45, and V49 (and a marginal detec- 
tion near V34). They are not detected in the shal- 
lower MG surveys and so no radio spectral index 
could be estimated. However, V35 has been de- 
tected in the Westerbor k Synthesis Radio Tele- 
scope (WSRT) survey of iTavlor et aD (|l996h . Its 
spectral index between 327 MHz and 1.4 GHz is 
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— 1.6, and is therefore very unlikely to be associ- 
ated with thermal emission. All of these sources 
are seen projected on Sh2-86, and so the radio 
"sources" could be just parts of the patchy ion- 
ization structure therein. As mentioned in SjSJ 
V39 is embedded at the end of a pillar adjacent 
to an ionization front caused by the exciting stars 
of NGC 6823, not V39. Inspection of the radio im- 
ages suggests that V45 and V49 might be similarly 
confused. 

Clearly, more sensitive radio observations, with 
spatial and spectral index information, are desir- 
able to further elucidate the true nature of all of 
these (luminous) sources. But a firm conclusion 
that should be emphasized is that none of these 
BLAST sources associated with NGC 6823 reveal 
a protostar or HMPO more massive than a BO to 
Bl star or about 20 Mq. There is thus no evi- 
dence for a new (or lagging) generation of stars 
that are as massive as the existing powerhouses of 
the NGC 6823 cluster, which reach about 40 M 



(07 V and 09.5 Ia; ISchaerer fc de Koterlll997h . 



7.3. Extended Submillimeter Emission 

In this paper we have focused on the compact 
sources in the BLAST05 maps of Vulpecula. More 
diffuse structure could also be st udied using al- 



gorit hms such as CLUMPFIND ([Williams et al 



1994). We leave such quantitative analysis for fu- 
ture work. However, it is worth illustrating a few 
features in the context of the discussion from ^5.21 
and the cool low luminosity-to-mass ratio clumps 
already identified. 

Between sources V02 and V03 there is a pair 
of extended BLAST features, neither sufficiently 
compact to be included in our source catalog. The 
first, to the west of the line, is a strong 100/im 
source, and even more prominent at 60/im; it is 
IRAS 19364+2252. This region exhibits a striking 
fan-shaped nebulosity in the IRAC 8 (im map (Fig- 
ure [17]). Since the object is also a radio source, it 
is probably an evolved object. It has weak emis- 
sion in the BLAST maps and is anti-coincident 
with 13 CO. The second source, right on the line, is 
possibly related to the IRA S source 19367+2251, 
which IGA/MIGA show to be much weaker as 
a compact source than the first. This extended 
BLAST source coincides with a peak in 13 CO and 
becomes brighter relative to the first source at in- 
creasing BLAST wavelengths. This behavior im- 



plies that it is a very cold object. 

Another example is associated with the pair of 
IRAS sources 19385+2245 and 19389+2233. They 
can be seen in all of the IRAS bands, but they 
are both diffuse in the BLAST maps. The first is 
weak in 13 CO and has complex structure in the 
8/im map. The second is elongated in BLAST 
and has a strong elongated CO cloud as well. It 
is also extended at 8 /xm, with a slight offset from 
the submillimeter source. 

A future analysis of such cold clouds in BLAST 
maps in conjunction with MIPS and IRAC source 
catalogs may be used to determine whether they 
represent an even earlier phase before condensa- 
tions form, or if they are simply distant cold en- 
velopes of established young stellar objects. 

8. CONCLUSIONS 

In this paper we present the first maps of the 
Galactic Plane observed with BLAST during its 
2005 LDB flight from Sweden to Canada, specif- 
ically the survey in the direction of star forming 
clouds in Vulpecula. 

Fourier deconvolution is used to improve the 
angular resolution of the maps from ~ 3f5 to 
~ 1' full-width half-power (the theoretical 500 /im 
diffraction limit). In these maps, 60 compact sub- 
millimeter sources are detected simultaneously at 
250, 350, and 500 /im. Complementary IRAS, 
MIPS GAL, and MSX photometry are used to con- 
strain the submillimeter-MIR SEDs of these ob- 
jects, and hence infer their cold dust temperatures 
and bolometric fluxes. Our sample has convinc- 
ingly revealed objects with a range of dust temper- 
atures, from ~ 12-40 K under the assumption of 
an isothermal modified blackbody with dust emis- 
sivity index /3 = 1.5. With these SED fits, we de- 
rive luminosity-to-mass ratios (independent of dis- 
tance) in the range 0.2-130 LqMq 1 . Those with 
low values are highly shielded and cool, with no 
evidence for star formation, while those with high 
values have embedded high-mass star formation. 

Using a 13 CO(l -> 0) data cube, the VLA 
Galactic Plane Survey H I cube, Spitzer IRAC 
and optical DSS images, we argue that 49 of the 
60 sources lie in a molecular cloud complex associ- 
ated with the open cluster NGC 6823 at ~ 2.3 kpc, 
10 objects are associated with the Perseus arm at 
~ 8.5 kpc and 1 object is in the outer Galaxy at 
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~ 14kpc. With these distance estimates we cal- 
culate bolometric luminosities and cloud masses 
associated with the thermal emission from cold 
dust. The most luminous object (7 x 10 4 L Q ) is 
in the Perseus arm. Near NGC 6823, the ranges 
are ~ 40-10 4 L Q , and ~ 15-700 M . 

A mass function is constructed for the 49 ob- 
jects associated with NGC 6823. It is compatible 
with the spectrum of molecular gas masses in other 
high-mass star forming regions, with a power-law 
index of -1.7 for the > 100 M Q sources. A flat- 
tening at lower masses might be present, but is 
affected by detection completeness of cool sources 
in this mass range. 

The luminosity-mass distribution we find is 
broadly consistent with an evolutionary sequence, 
from cool high-mass, low-luminosity clumps (most 
are not detected with IRAS or MSX and there 
is no evidence for any star formation), to more 
evolved HMPOs and UC H II regions. But in- 
terestingly, among the embedded objects in this 
molecular complex near NGC 6823 - the next gen- 
eration - there are none quite as massive as the 
40 M Q stars currently powering the nebula. 
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Table 2 
BLAST Sources 



BLAST 


Source name 


^250 


0250 


^350 


C350 


-F500 


C500 


Vlsr 


ID 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


km s _1 


V01 


BLAST 


J193837+230541 


59.9 


4.2 


27.1 


2.8 


11.0 


1.2 


30.8 


V02 


BLAST 


J193849+230839 


258.0 


5.1 


97.8 


3.4 


36.9 


1.3 


32.3 


V03 


BLAST 


J193858+224637 


197.9 


4.7 


101.9 


3.5 


44.5 


1.3 


36.5 


V04 a 


BLAST 


J193921+224401 


23.6 


4.1 


8.9 


2.7 


7.3 


1.2 


27.3 


V05 


BLAST 


J193935+235947 


325.7 


5.6 


124.1 


3.8 


43.6 


1.3 


37.0 


V06 


BLAST 


J193943+232601 


20.4 


4.1 


13.1 


2.7 


5.8 


1.2 


43.2 


V07 b 


BLAST 


J194046+241007 


38.7 


4.1 


12.2 


2.7 


5.2 


1.1 


-53.5 


V08 


BLAST 


J194059+240439 


225.6 


4.9 


81.6 


3.2 


30.7 


1.2 


34.8 


V09 


BLAST 


J194103+240111 


194.8 


4.7 


102.1 


3.5 


41.2 


1.3 


34.4 


V10 


BLAST 


J194106+235513 


21.5 


4.1 


12.6 


2.7 


6.7 


1.2 


34.0 


Vll 


BLAST 


J194136+232325 


18.3 


4.0 


18.9 


2.7 


10.1 


1.2 


27.8 


V12 


BLAST 


J194144+231410 


33.2 


4.1 


16.9 


2.7 


7.1 


1.1 


29.4 


V13 


BLAST 


J194144+231804 


26.6 


4.1 


16.4 


2.7 


7.5 


1.2 


23.7 


V14 


BLAST 


J194146+234914 


19.0 


4.1 


10.3 


2.7 


3.3 


1.1 


34.6 


V15 


BLAST 


J194213+232305 


41.9 


4.1 


21.8 


2.7 


10.5 


1.2 


24.6 


V16 


BLAST 


J194217+231958 


23.5 


4.1 


10.8 


2.7 


5.5 


1.2 


33.0 


VI T 


BLAST 


J194226+234532 


30.4 


4.1 


10.8 


2.7 


5.6 


1.2 


-2.1 


V18 


BLAST 


J194228+230458 


259.9 


5.1 


114.3 


3.6 


42.8 


1.3 


26.6 


V19 ac 


BLAST 


J194232+234913 


15.8 


4.1 


5.6 


2.7 


2.2 


1.1 


-3.0 


V20 c 


BLAST 


J194232+234253 


27.6 


4.1 


17.9 


2.7 


9.3 


1.2 


0.8 


V21 c 


BLAST 


J194234+234412 


20.9 


4.1 


15.6 


2.7 


6.6 


1.2 


-0.3 


V22 


BLAST 


J194234+232321 


60.9 


4.1 


34.8 


2.8 


15.1 


1.2 


30.5 


V23 c 


BLAST 


J194236+235003 


129.6 


4.4 


45.8 


2.9 


15.7 


1.2 


-3.0 


V24 c 


BLAST 


J194244+234021 


62.5 


4.2 


37.1 


2.8 


12.3 


1.2 


0.3 


V25 


BLAST 


J194246+225536 


30.1 


4.1 


17.7 


2.7 


10.3 


1.2 


26.9 


V26 


BLAST 


J194254+232408 


73.2 


4.2 


42.4 


2.8 


18.8 


1.2 


31.4 


V27 


BLAST 


J194306+230125 


19.8 


4.1 


14.1 


2.7 


7.0 


1.2 


27.0 


V28 


BLAST 


J194308+232457 


21.0 


4.1 


15.0 


2.7 


6.5 


1.2 


33.1 


V29 a 


BLAST 


J194311+232010 


15.7 


4.1 


2.3 


2.7 


5.5 


1.2 


24.0 


V30 


BLAST 


J194311+234405 


852.7 


11.0 


378.5 


8.5 


153.4 


2.7 


22.8 


V31 


BLAST 


J194311+232237 


13.3 


4.1 


19.1 


2.7 


2.3 


1.1 


33.4 


V32 


BLAST 


J194317+231352 


221.2 


4.9 


110.7 


3.6 


43.3 


1.3 


29.3 


V33 


BLAST 


J194319+232639 


17.5 


4.0 


17.7 


2.7 


6.7 


1.2 


22.5 


V34 a 


BLAST 


J194319+233906 


27.5 


4.1 


24.1 


2.7 


7.4 


1.2 


22.9 


V35 


BLAST 


J194323+231512 


22.0 


4.1 


7.9 


2.7 


1.4 


1.1 


30.4 


V36 a 


BLAST 


J194325+234503 


86.6 


4.2 


68.3 


3.1 


18.9 


1.2 


21.3 


V37 


BLAST 


J194326+240618 


24.8 


4.1 


9.3 


2.7 


3.5 


1.1 


20.7 


V38 a 


BLAST 


J194327+235638 


64.4 


4.2 


43.8 


2.9 


22.6 


1.2 


22.4 d 


V39 


BLAST 


J194328+232032 


57.8 


4.1 


31.1 


2.8 


12.6 


1.2 


25.7 


V40 


BLAST 


J194329+234013 


189.7 


4.7 


86.9 


3.3 


35.4 


1.3 


20.5 


V41 


BLAST 


J194331+235728 


90.8 


4.2 


42.7 


2.8 


16.3 


1.2 


22.8 


V42 a 


BLAST 


J194333+231445 


23.8 


4.1 


24.0 


2.7 


7.5 


1.2 


26.9 


V43 a 


BLAST 


J194334+232737 


41.8 


4.1 


24.7 


2.8 


3.6 


1.1 


26.9 


V44 


BLAST 


J194335+235332 


28.2 


4.1 


14.7 


2.7 


5.5 


1.2 


21.4 
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This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Table 2 — Continued 



BLAST 


Source name 


i*250 


0250 


-F350 


0350 


F500 


0500 


V\ at 


ID 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


km s _1 


V45 a 


BLAST 


J194338+231725 


10.6 


4.1 


7.0 


2.7 


0.4 


1 


.1 


26.6 


V46 


BLAST 


J194340+235506 


28.8 


4.1 


12.5 


2.7 


5.3 


1 


.1 


22.5 


V47 


BLAST 


J194342+232006 


129.9 


4.4 


70.2 


3.1 


33.8 


1 


.3 


26.7 


V48 


BLAST 


J194344+240005 


23.5 


4.1 


11.2 


2.7 


3.3 


1 


.1 


22.2 


V49 a 


BLAST 


J194346+232532 


27.0 


4.1 


17.9 


2.7 


11.5 


1 


.2 


26.1 


V50 


BLAST 


J194347+233505 


18.6 


4.1 


10.1 


2.7 


3.0 


1 


.1 


22.1 


V51 c 


BLAST 


J194349+230840 


22.4 


4.1 


16.1 


2.7 


6.8 


1 


.2 


-7.0 


V52 


BLAST 


J194350+232839 


486.4 


7.1 


214.5 


5.3 


84.2 


1 


.8 


27.6 


V53 


BLAST 


J194350+231846 


40.6 


4.1 


20.0 


2.7 


9.1 


1 


.2 


27.7 


V54 c 


BLAST 


J194351+234625 


21.8 


4.1 


7.0 


2.7 


2.4 


1 


.1 


3.2 C 


V55 a 


BLAST 


J194402+232409 


30.6 


4.1 


21.9 


2.7 


5.0 


1 


.2 


29.3 


V56 


BLAST 


J194404+234418 


53.5 


4.1 


37.3 


2.8 


15.6 


1 


.2 


21.7 


V57 


BLAST 


J194405+232658 


40.7 


4.1 


19.6 


2.7 


7.1 


1 


.2 


27.2 


V58 c 


BLAST 


J194409+230437 


15.4 


4.1 


10.7 


2.7 


5.2 


1 


.2 


-12.7 


V59 c 


BLAST 


J194424+234239 


23.8 


4.1 


12.9 


2.7 


4.1 


1 


.1 


-10.7 


V60 


BLAST 


J194430+233426 


14.9 


4.1 


14.1 


2.7 


4.8 


1 


.2 


26.2 



Note. — Flux densities for BLAST sources are quoted at precisely 250, 350, and 500 /zm us- 
ing SED fits to obtain color-corrections for the band-averaged flux densities (<2]). The quoted 
statistical uncertainties are determined from Monte Carlo simulations (£ )2.6p , and do not include 
calibration uncertainties. 

a These sources are located on ripples in the deconvolved map and the BLAST colors are con- 
sidered unreliable as a result. 

b V07 is believed to lie in the outer galaxy. 

c These sources are associated with a molecular cloud in the Perseus arm. 
d Also has a comparable component at — 15kms _1 . 
c Also has a comparable component at 23kms~ 1 . 
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Table 3 
IRAS PSC Counterparts 



BLAST 


IRAS 


Aa 


a 


AS 1 




F%2 


F25 


-^60 


-F100 




ID 


ID 


(") 




(") 




(Jy) 


(Jy) 


(Jy) 


(Jy) 




V01 


19364+2258 


-19. 


.4 


-3 


.3 


0.8 


8.3 


389.0 b 


698 


,0 b 


V02 


19366+2301 


-21 


.0 


-1. 


.0 


8.8 


62.0 


389.0 


698 


.0 


V03 


19368+2239 


-4 


.5 


-5. 


.4 


2.3 


11.7 


112.0 


231 


.0 


V05 


19374+2352 


-27 


.0 


7. 


.3 


4.6 


24.3 


422.0 


768 


.0 


V07 


19386+2403 


-5 


.9 


2. 


.1 


1.8 


4.7 


71.2 


63 


.8 


V08 


19388+2357 


-3 


.0 


-0. 


.1 


1.3 


10.9 


254.0 


433 


.0 


V09 


19389+2354 


1. 


.5 


-0. 


.1 


1.8 


1.2 b 


58.2 


433 


,0 b 


V14 


19396+2342 


3. 


.1 


9 


.6 


1.5 b 


0.8 b 


3.1 


16. 


.7 


V18 


19403+2258 


-21 


.0 


13. 


.9 


16.2 


118.0 


562.0 


640 


.0 


V90 

V iU 




ft 
— U 


n 


n 


9 


1 Q b 


1 (i 

1 . u 


A Q b 


00. 


. U 


V23 c 


19404+2342 


-24 


.0 


-7. 


.9 


9.7 


35.3 


439.0 


980 


.0 


V24 


19406+2333 


0. 


.0 


8. 


.1 


0.7 


5.1 


35.0 


89. 


.8 


V26 


19407+2316 


-4 


.5 


-6. 


.8 


4.1 b 


4.2 


23.2 


165 


,0 b 


V30 


19410+2336 


4. 


.5 


0. 


.9 


14.4 


109.0 


982.0 


1630 


.0 


V32 


19411+2306 


12 


.0 


6. 


.2 


5.3 


41.2 


159.0 


557 


.0 


V37 


19413+2358 


-3 


.0 


-8. 


.7 


15.4 b 


0.8 b 


4.5 


37 


,0 b 


V38 d 


19413+2349 


10 


.5 


20. 


.1 


0.7 


4.6 


42.5 


132 


.0 


V40 


19413+2332 


-4 


.5 


-9. 


.3 


4.2 


38.5 


217.0 


1630 


,0 b 


V51 


19416+2301 


0. 


.0 


-12 


.9 


1.0 b 


1.0 


4.7 b 


35. 


,8 b 


V52 


19416+2321 


-18 


.0 


20 


.5 


5.0 b 


24. 7 b 


268.0 


624 


.0 


V53 d 


19416+2312 


-39 


.0 


39. 


.3 


0.6 


10.7 b 


13.6 b 


250 


.0 


V57 


19419+2319 


4. 


.5 


-15 


.0 


2.8 


2.6 


36. l b 


624 


,0 b 


V59 


19422+2335 


-1 


.5 


2 


.1 


1.4 


2.3 


25.6 


46. 


.6 



Note. — IRAS point sources associated with BLAST objects. The search 
radius is a variable function of both BLAST and IRAS positional uncertain- 
ties (see £)3.ip . 

a BLAST source tangent plane offsets (E and N) compared to IRAS source 
position. 

b These IRA S measurements are flagged as upper-limits in the catalog (usu- 
ally due to confusion). Values from Table H] are used instead. 

C IRAS 19404+2342 is close to both V19 and V23; the ambiguity is resolved 
with the improved resolution of the MIPS 70 /im map. 

d Although IRAS 19413+2349 and 19416+2312 are close to V38 and V53 
respectively, their SEDs are not compatible with the BLAST photometry 
and the identifications are not used in subsequent analyses. 
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Table 4 

IRAS and Spitzer MIPS Photometry 



BLAST 


F 12 


Flag 


F 25 


Flag 


-^60 


Flag 


-Fioo 


Flag 


F 70 


Flag 


ID 


(Jy) 




(Jy) 




(Jy) 




(Jy) 




(Jy) 




V01 


1.2 a 


g 


7.2 a 


g 


46.5 


g 


120.5 


u 


33.0 


s 


V02 


124 a 


g 


60.9 a 


g 


437.7 a 


g 


733. 6 a 


g 


211.1 


s 


V03 


2.4 a 


g 


12.4 a 


g 


114.5 a 


g 


243. 2 a 


g 


78.2 


s 


V04 


0.2 


n 


0.2 


n 


0.2 


n 


31.3 


n 


0.4 


n 


V05 


5.9 a 


g 


28.6 a 


g 


410. 5 a 


g 


772. 9 a 


g 


218.6 


s 


V06 


-0.1 


n 


-0.1 


11 


4.2 


g 


24.7 


u 


2.5 


g 


V07 


1.7 a 


g 


4.7 a 


g 


66.7 a 


g 


135. 7 a 


g 


50.4 


s 


V08 


2.3 a 


g 


12. l a 


g 


259. 9 a 


g 


569. 6 a 


g 


206.7 


s 


V09 


1.5 a 


g 


2.9 


g 


44.1 a 


g 


256.8 


u 


44.1 


s 


V10 


0.0 


n 


0.0 


g 


2.4 


11 


48.3 


u 


0.9 


n 


Vll 


-0.1 


n 


-0.2 


11 


-0.2 


n 


18.2 


11 


0.8 


n 


V12 


1.7 


g 


4.0 


g 


14.8 


g 


70.9 


u 


9.4 


g 


V13 


-0.4 


n 


-0.4 


11 


-3.6 


n 


13.0 


11 


-1.5 


n 


V14 


0.3 


n 


0.5 


u 


4.0 a 


g 


28.0 a 


u 


2.7 


g 


V15 


-0.2 


n 


0.2 


11 


9.4 


11 


83.3 


u 


3.1 


n 


V16 


1.2 


u 


2.0 


n 


11.7 


u 


91.9 


u 


3.5 


g 


V17 


0.2 


n 


0.9 


u 


5.9 


u 


14.5 


11 


3.4 


n 


V18 


17.2 a 


g 


108. 5 a 


g 


570. a 


g 


699. l a 


g 


275.4 


s 


V19 


17.6 


u 


52.8 


u 


460.1 


11 


979.5 


u 


148.6 


s 


V20 


0.0 


g 


0.3 a 


g 


6.6 


11 


19.4 


11 


6.3 


g 


V21 


-0.1 


n 


-0.3 


11 


-1.9 


n 


-0.1 


n 


4.5 


g 


V22 


0.7 


n 


0.5 


11 


7.4 


11 


97.7 


u 


6.0 


n 


V23 


15.3 a 


u 


53.3 a 


u 


506.5 a 


g 


916.8 a 


g 


192.2 


s 


V24 


1.4 a 


u 


3.5 a 


g 


32.0 a 


g 


110. 6 a 


u 


25.5 


s 


V25 


0.2 


n 


0.9 


u 


3.8 


g 


47.0 


g 


0.7 


g 


V26 


2.0 


u 


4.0 a 


g 


21.9 a 


g 


91.4 


u 


7.7 


g 


V27 


0.4 


u 


0.7 


u 


7.7 


u 


51.7 


u 


0.7 


n 


V28 


2.6 


u 


4.2 


u 


-3.3 


n 


36.9 


11 


3.2 


n 


V29 


-0.7 


n 


0.1 


11 


-0.5 


n 


-6.2 


n 


-0.9 


g 


V30 


18.5 a 


g 


118.0 a 


g 


892. 2 a 


g 


1463.9 a 


g 


378.2 


s 


V31 


0.7 


n 


0.9 


u 


2.0 


n 


-7.8 


n 


-0.3 


n 


V32 


6.1 a 


g 


37.3 a 


g 


151. 2 a 


g 


478. 2 a 


g 


62.6 


s 


V33 


1.2 


n 


1.9 


u 


21.5 


u 


37.0 


u 


-1.3 


n 


V34 


3.2 


u 


5.4 


u 


83.9 


u 


334.1 


u 


2.2 


n 


V35 


6.2 


u 


22.9 


u 


87.9 


u 


461.9 


u 


6.6 


n 


V36 


5.6 


u 


8.6 


g 


88.4 


11 


472.2 


u 


42.6 


s 


V37 


0.3 


n 


0.8 


n 


4.4a 


g 


24.6 


u 


3.3 


g 


V38 


0.9 


u 


5.3 


u 


34.0 


11 


120.9 


u 


22.9 


s 


V39 


3.8 


u 


3.7 


g 


40.4 


11 


183.2 


u 


19.4 


g 


V40 


7.2 a 


g 


44.4a 


g 


214.3 a 


g 


400.5 


g 


152.3 


s 


V41 


0.6 


n 


5.2 


u 


33.5 


u 


131.7 


u 


17.8 


u 


V42 


-0.8 


n 


-0.6 


n 


-0.6 


n 


223.5 


u 


-3.3 


n 


V43 


1.8 


n 


1.6 


11 


21.1 


u 


147.6 


u 


7.8 


n 


V44 


0.0 


n 


0.1 


11 


3.6 


11 


35.2 


n 


2.0 


u 
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Table 4 — Continued 



BLAST 


F\2 


Flag 


-Fas 


Flag 


-^60 


Flag 


-F100 


Flag 


-F70 


Flag 


ID 


(Jy) 




(Jy) 




(Jy) 




(Jy) 




(Jy) 




V45 


0.0 


n 


3.9 


11 


1.7 


11 


17.7 


11 


6.8 


n 


V46 


0.0 


n 


-0.2 


11 


9.5 


u 


137.3 


11 


12.3 


11 


V47 


0.9 


g 


10.5 


g 


60.5 


g 


174.6 


g 


21.0 


s 


V48 


0.8 


n 


1.2 


11 


9.0 


u 


61.0 


11 


8.2 


g 


V49 


-0.9 


n 


1.0 


11 


16.7 


u 


149.0 


11 


17.6 


g 


V50 


0.2 


g 


1.0 


g 


6.7 


u 


51.5 


11 


4.2 


g 


V51 


0.3 


n 


1.9 a 


11 


15.9 


11 


103.5 


11 


3.2 


g 


V52 


4.4 


u 


27.6 


11 


278. 2 a 


g 


675. l a 


g 


211.7 


s 


V53 


-0.5 a 


n 


-0.3 


11 


5.5 


11 


227.9 a 


11 


0.7 


n 


V54 


0.5 


u 


0.9 


u 


11.2 


11 


82.4 


11 


12.9 


g 


V55 


1.9 


u 


5.6 


11 


11.3 


g 


369.6 


11 


8.1 


n 


V56 


-0.1 


n 


-0.3 


11 


1.1 


11 


31.4 


11 


-0.4 


n 


V57 


1.7 a 


g 


5.3 a 


g 


4.5 


11 


423.2 


11 


5.7 


g 


V58 


-0.1 


n 


-0.2 


11 


7.2 


11 


90.1 


11 


2.5 


g 


V59 


1.4 a 


g 


2.3 a 


g 


26.4 a 


g 


57.4 a 


g 


24.7 


s 


V60 


0.8 


n 


1.3 


u 


10.6 


11 


75.4 


11 


2.2 


n 



Note. — IRAS and Spitzer MIPS map photometry are described in §3. H and i )3.2l respec- 
tively. Each measurement is flagged upon visual inspection: "g" if it is good; "u" for upper 
limit; "n" for no visible source; and "s" for saturated (unreliable) values in the MIPS 70 ^m 
map. 

a Clear detections from the IRAS PSC (Table [3]) are used in favor of these values when 
available. 
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Table 5 
msx counterparts 



BLAST 


Aa a 


A5 a 


F s .3 


^12.1 


Fl4.7 


-F21.3 


ID 


(") 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


V01 


-6.9 


1.4 


0.9 


0.8 


0.8 


4.7 


V02 


-11.0 


-1.3 


4.3 


5.3 


3.7 


23.9 


V03 


18.0 


8.1 


0.9 


1.1 


0.6 


1.6 




-20.7 


-9.9 


0.8 


1.0 


1.4 


3.4 


V05 


-6.9 


26.0 


1.7 


1.8 


1.2 


2.7 




-2.7 


2.9 


0.9 


0.8 


1.2 


6.2 


V06 


19.3 


-63.1 


0.4 








V07 


-0.0 


11.2 


1.0 


1.4 


0.8 


2.0 




-26.0 


-53.6 


0.1 


1.2 






V08 


-0.0 


-23.9 


0.6 


1.4 


0.9 






-8.2 


6.7 


0.5 


1.3 


0.8 


3.6 


V12 


-48.2 


-41.7 


1.7 


1.2 


0.7 






8.3 


9.8 


0.5 


2.0 


3.1 


4.7 


V13 


-4.1 


21.1 


0.1 








V14 


-11.0 


16.7 


0.2 




0.8 




V18 


-8.3 


16.9 


8.4 


9.0 


6.2 


57.5 


V19 


11.0 


28.5 


0.9 


1.8 


2.2 


3.9 




37.0 


38.2 


0.7 


2.6 


3.2 


7.9 




-9.6 


49.0 


0.1 


0.6 




1.3 


V20 


-0.0 


-0.2 


0.2 








V21 


-15.1 


13.9 


0.1 








V22 


-20.7 


-3.6 


0.2 








V23 


-15.1 


-11.7 


0.7 


2.6 


3.2 


7.9 




17.8 


-5.2 


1.5 


2.5 


1.3 


3.1 




-19.2 


17.1 


0.8 


2.2 


2.7 


5.2 


V24 


-5.5 


9.9 


0.4 






2.5 


V26 


15.1 


6.7 


1.4 


2.1 


2.2 


3.6 


V28 


-23.4 


-28.5 


0.9 


1.4 


1.7 


1.9 


V30 


-0.0 


-0.2 


5.3 


8.5 


14.9 


47.4 


V31 


-15.1 


17.8 


0.2 








V32 


8.3 


8.5 


1.8 


5.2 


8.7 


25.6 


V35 


66.1 


25.8 


0.2 






2.7 


V36 


31.6 


-15.4 


1.7 


1.6 


1.0 


2.1 


V37 


20.5 


60.9 


0.9 


0.6 


0.5 




V38 


2.7 


16.0 


0.7 


0.9 


1.2 


2.7 


V39 


2.8 


10.5 


0.9 


1.9 


2.1 


2.3 


V40 


-9.6 


-21.8 


2.2 


2.5 


1.1 


6.1 




6.9 


11.7 


2.0 


2.2 


2.7 


18.3 


V43 


-6.9 


-20.0 


0.2 








V46 


8.2 


4.9 


0.1 










4.1 


-39.0 


0.3 


1.8 






V47 


4.1 


8.8 


0.1 


1.1 


1.7 


4.7 


V49 


-19.3 


-24.1 


0.4 


0.8 






V50 


15.1 


48.0 


0.2 









Table 5 — Continued 



BLAST 


Aa a 


AS a 


F 8 .3 


-F12.1 


-F14.7 


-P21.3 


ID 


(") 


(") 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


V51 


-2.8 


-7.9 


0.2 








V52 


-4.1 


0.0 


0.2 






3.3 


V53 


-20.7 


-29.8 


0.1 








V55 


-17.9 


11.8 


0.2 








V57 


-1.4 


-11.7 


1.3 


1.3 


1.2 


1.2 


V59 


1.4 


1.5 


0.8 


1.1 


0.6 




V60 


-52.2 


5.7 


0.1 


1.0 







Note. — All candidate MSX counterparts are listed for each 
BLAST source (§ 



a BLAST source tangent plane offsets (E and N) compared 
to MSX source position. 
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Table 6 
Results from SED Fits 



BLAST 


Temperature 






Luminosity Fraction 


L F m/M c 
(LqMq 1 ) 


ID 


(K) 


(M ) 


(L©) 


(%) 


V01 


27.2 ± 1.5 


44 ±7 


750 ± 130 


76 


17.5 ±5.0 


V02 


31.9 ±2.3 


114 ±22 


4500 ± 1000 


78 


40.5 ± 15.8 


V03 


22.2 ±0.9 


273 ± 40 


1470 ± 150 


67 


5.6 ± 1.2 


V04 


20.0 


38 ±5 


121 ±15 




3.2 


V05 


31.3 ± 1.9 


130 ±27 


4910 ±880 


89 


36.1 ±12.2 


V06 


20.2 ±2.4 


55 ±35 


96 ±57 


62 


1.9 ± 1.9 


V07 


32.6 ± 1.9 


579 ± 122 


28200 ±4900 


83 


49.0 ± 16.5 


V08 


29.7 ± 1.7 


97 ± 19 


2610 ±430 


84 


26.5 ±8.8 


V09 


22.8 ± 1.7 


233 ± 45 


1560 ±350 


80 


6.4 ±2.4 


V10 


15.8 ±2.6 


89 ±50 


53 ±22 


93 


0.9 ±0.9 


Vll 


12.3 ± 1.6 


213 ± 125 


42 ± 15 




0.2 ±0.1 


V12 


21.9 ± 1.0 


42 ±8 


220 ± 30 


51 


5.3 ± 1.3 


V13 


16.5 ± 1.5 


85 ±28 


87 ±20 


70 


1.0 ±0.5 


V14 


22.6 ± 1.0 


17 ±5 


107 ± 11 


70 


6.0 ± 1.6 


V15 


17.8 ± 1.5 


93 ±23 


145 ± 43 




1.5 ±0.7 


V16 


20.2 ±2.0 


39 ± 16 


117 ±34 




2.9 ± 1.6 


V17 


19.0 ± 1.1 


651 ± 171 


1400 ± 280 




2.3 ±0.8 


V18 


30.3 ± 1.5 


148 ±24 


4590 ± 740 


58 


30.9 ±9.3 


V19 


20.0 


281 ±67 


890 ±210 




3.2 


V20 


16.6 ±3.1 


1274 ± 739 


1250 ±590 


75 


0.9 ±0.9 


V21 


15.8 ±3.2 


1109 ± 762 


920 ± 540 


66 


0.9 ±0.9 


V22 


17.7 ± 1.3 


145 ± 34 


227 ± 44 


72 


1.6 ±0.6 


V23 


39.4 ± 1.7 


526 ± 65 


71000 ± 11000 


89 


130.4 ±31.6 


V24 


23.2 ±0.8 


970 ± 121 


6770 ± 640 


68 


6.7 ± 1.0 


V25 


17.6 ±2.6 


93 ±44 


112 ±61 


100 


1.3 ± 1.1 


V26 


18.9 ± 1.1 


150 ±27 


346 ± 62 


58 


2.4±0.7 


V27 


14.9 ±2.6 


105 ±64 


54 ±26 




0.6 ±0.6 


V28 


15.9 ± 1.8 


90 ±46 


60 ±25 


39 


0.8 ±0.7 


V29 


20.0 


23 ±4 


73 ± 14 




3.2 


V30 


25.9 ± 1.0 


680 ± 115 


9020 ± 700 


63 


13.0 ±2.9 


V31 


16.8 ± 1.7 


43 ± 14 


45 ±22 


50 


1.0 ±0.6 


V32 


26.2 ± 1.3 


208 ± 30 


2820 ±410 


69 


13.7 ±3.7 


V33 


14.2 ±2.3 


107 ±64 


49 ±21 




0.4 ±0.4 


V34 


20.0 


52 ±6 


166 ± 18 




3.2 


V35 


20.7 ± 1.2 


22 ±7 


84 ± 19 


54 


4.0 ± 1.1 


V36 


20.0 


118 ± 12 


374 ± 39 


54 


3.2 


V37 


22.6 ±0.9 


24 ±6 


139 ± 18 


73 


5.8 ± 1.5 


V38 


20.0 


132 ± 13 


417 ±42 


72 


3.2 


V39 


21.7±0.6 


78 ±8 


390 ± 32 


70 


5.0 ±0.7 


V40 


28.0 ± 1.6 


142 ± 26 


2670 ±610 


66 


18.3 ±6.9 


V41 


20.8 ±0.8 


110 ± 14 


423 ± 43 




3.8 ±0.8 


V42 


20.0 


49 ±6 


155 ± 18 




3.2 


V43 


20.0 


45 ±5 


142 ± 15 


86 


3.2 


V44 


18.3 ± 1.3 


53 ± 15 


103 ±21 




2.0 ±0.7 
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Table 6 — Continued 



BLAST 


TpTrmPTfttnrr 1 
XClAlLlCx at 111 c 


M 


^FIR 


1 ,il TT11 nAQl i~\T HrflP'fin'tl 
J_J LlllllllUol Li y _L 1 CtHjlUll 


-^FIR/ lvl c 

(LqMq 1 ) 


ID 


(K) 


(Ms) 


(Lo) 


(%) 


V45 


20.0 


14 ±4 


46 ± 13 




3.2 


V46 


21.4 ±1.7 


37 ± 12 


163 ±35 


66 


4.5 ± 1.4 


V47 


20.2 ±2.9 


219 ±70 


830 ±310 


61 


3.4 ±2.3 


V48 


22.5 ±1.7 


25 ±7 


154 ± 34 




5.8 ±2.3 


V49 


20.0 


64 ±6 


203 ± 18 


58 


3.2 


V50 


21.1 ±2.1 


25 ± 10 


101 ±31 


72 


4.0 ±2.4 


V51 


16.6 ±3.4 


1088 ± 727 


1050 ±570 


70 


1.2 ± 1.2 


V52 


24.6 ± 1.4 


390 ± 75 


3850 ±610 


84 


9.7±3.1 


V53 


18.5 ±1.3 


80 ± 17 


156 ± 34 


80 


1.8 ±0.7 


V54 


24.9 ±1.3 


210 ±48 


2320 ±410 




10.5 ±3.2 


V55 


20.0 


48 ±6 


153 ±17 


59 


3.2 


V56 


16.1 ±1.3 


180 ± 42 


176 ± 49 




0.9 ±0.5 


V57 


20.1 ±1.1 


58 ±9 


181 ±31 


60 


3.1 ± 1.0 


V58 


16.7 ±3.9 


757 ±510 


540 ± 330 




0.8 ±0.7 


V59 


25.4 ±1.5 


297 ±67 


3400 ± 390 


63 


12.2 ±4.4 


V60 


15.8 ±2.3 


65 ±28 


42 ±21 


34 


0.6 ±0.5 



Note. — Results from modified blackbody SED (Equation 2} fits to submillimcter and 
FIR photometry (j]IJ). These fits assume a dust emissivity index (3 = 1.5. Cloud masses M c 
are derived from the amplitude A of the SED using Equation [5] and distances from <J5] FIR 
luminosities, Lfir, are the integrated luminosities from the modified blackbody. Objects 
with MIR photometry have estimates for the total bolomctric luminosity, which may be 
obtained by dividing the FIR luminosity by the Luminosity Fraction. Quoted uncertainties 
are statistical errors derived from Monte Carlo simulations f£ )4.2p . Other systematic errors are 
not included. For example, adopting /3 = 2.0 decreases temperatures by ~ 5K, and increases 
cloud masses by a factor ~ 2 (JOJ. There are also significant distance uncertainties (fj5]), 
and a factor of ~ 2 uncertainty introduced through choices of kq and R in Equation [5] 
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